BROWNIAN MOTIONS ON METRIC GRAPHS: 
FELLER BROWNIAN MOTIONS ON INTERVALS REVISITED 
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Dedicated to the memory of Pierre Duclos 

ABSTRACT. The construction of the paths of all possible Brownian motions (in the 
sense of |21| ) on a half line or a finite interval is reviewed. 
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1. Introduction 

In recent years there has been a growing interest in metric graphs as the underlying 
structure for models in many branches of science, for example in physics, biology, 
chemistry, engineering and computer science, to name just a few. The interested 
reader is referred to the review [31 1 and to the articles in the volume [4]. Metric graphs 
are piecewise linear varieties with singularities at a finite number of points, namely at 
the vertices, and they can be thought of as a finite collection of finite intervals or half 
lines, which are glued together at some of their endpoints. 

There exists an extensive amount of literature on Laplace operators and their semi- 
groups on metric graphs, see e.g., [|22][23][27]-[29] and the literature quoted there. In 



view of this, it is natural to investigate Brownian motions on metric graphs, and in this 
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context we also want to mention the articles [2, 5j [T3] - [T6l[30l which deal with various 
aspects of stochastic processes on (metric) graphs. 

The present article is one of four articles (cf. (24]-[26j) of the authors in which 
we address the problem of the characterization of all Brownian motions (the precise 
definition of this class of stochastic processes is given in [24J), and their pathwise 
construction on a given metric graph. Thus, for metric graphs we consider — with 



certain restrictions (see |24| and below) — the analogue of the problem raised by 



Feller in his pioneering articles |9|4T2[, which have to be considered together with the 
work of Ito and McKean ]T8|[T9| 7 and the article [38 ] of Wentzell. For later accounts 
of this subject we also refer to J7J and (2TJ. 

Of course, according to our above description, the simplest metric graphs are given 
by a half line as M + or a compact interval, say, [0, 1], that is, the cases considered in 
the aforementioned classical literature. One of the intentions of the present article is 
to bring the material presented there into a form which is suitable for a generalization 
to metric graphs. The other intention is to give a pedagogical (and rather detailed) 
account of the construction of Feller Brownian motions on intervals. A very readable 
treatment can be found in the book [21 ] by Knight. However, there some arguments 
are only hinted at, others are not so easy to follow (at least for the present authors). 
We have tried to make the article self-contained, which means that several of our argu- 
ments are rather well-known. On the other hand, we shall give a number of arguments 
and computations which to the best of our knowledge are new. For example, we shall 
compute all transition and resolvent kernels explicitly, and in two cases we provide 
rather simple calculations of the generators based on Dynkin's formula (61. These 
calculations carry through to the case of metric graphs, and this will make it possible 



to compare the resolvents obtained with those found in |27 28]. This in turn allows 
— at least in certain cases — to rediscover from a stochastic point of view a central 
entity of quantum theory, namely the scattering matrix. 

For the definition of a Brownian motion on an interval we shall follow Knight (2TJ 



(in 1 24 1 this definition is extended to metric graphs). Let I be a finite or semi-finite 
interval, and let X denote a stochastic process with values in / U {A}, where A is 
a cemetery point. Moreover, let B be a standard one dimensional Brownian motion. 
We denote by X, B, the processes X, B respectively, with absorption (i.e., stopping) 
at the endpoint(s) of I. 

Definition 1.1. A Brownian motion X on a finite or semi-infinite interval is a normal 
strong Markov process on IU {A}, a.s. with right continuous paths, continuous paths 
up to its lifetime, and such that the stochastic process X is equivalent to B. 

We remark that this definition implies that we consider a considerably smaller class 
of stochastic processes than has been treated in the quoted work by Feller and Ito- 
McKean, in that there the paths are not required to be continuous up to the lifetime. 
In particular, there jumps from an endpoint back into the interval are allowed. On 
the other hand, this restriction will allow us in the present series of articles to work 
within the class of Feller processes, which has advantages concerning the control of 
the strong Markov property. Some of the consequences the restriction we impose will 
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also be discussed below. Within the more general framework of metric graphs the 
assumption that the paths are continuous up to their lifetime will be removed in a 
forthcoming work. 

For the sake of definiteness, from now on we shall only consider intervals I of 
the form [0, +00) or [0, 1]. Co([0, +00)) denotes the space of real valued continuous 
functions on [0, +00) vanishing at infinity, while C([0, 1]) is the space of real valued 
continuous functions on [0, 1]. By Co (J) we mean either of the two spaces, and 
we make the convention that every / G Co (J) is extended to I U {A} by setting 
/(A) = 0. Cq(I) denotes the subspace of Cq(I) consisting of those functions / in 
Co (I) which are twice continuously differentiable in the interior of /, and are such 
that the second derivative /" extends to a continuous function on I. It is not difficult 
to check that for / G C5 2 (I) the derivative /' has a (finite) limit /(0+) from the right 
at 0, and — if I = [0, 1] — also a (finite) limit f(l— ) from the left at 1. Then for 
Brownian motions on / as in definition | 1 . 1 | Feller 's theorem is the following statement. 

Theorem 1.2. Assume that X is a Brownian motion on I. Then the generator A of 
its semigroup on Cq(I) is given by Af = 1/2/" with a domain D(A) being a subset 
o/Cq(J). Moreover the following holds true: 

(a) Suppose that I = [0, +00). Then there exist ao, bo, Co G [0, 1] with ao + 60 + 
Co = 1, Qo 7^ 0, such that every f G T>(A) satisfies the Wentzell boundary 
condition at the origin 

(1.1) ao/(0) - bof(0+) + I /"(0+) = 0, 

and the domain T>(A) C Cq([0, +00)) is uniquely determined by this bound- 
ary condition. 

(b) Suppose that I = [0, 1]. Then there exist a v , b v , c v G [0, 1], v = 0, 1, with 
a v + b v + c v = 1, a v 7^ 0, such that every f G T>(A) satisfies the boundary 
condition ( |1.1[ ), supplemented by the boundary condition 

(1-2) ai/(l) + 6i/'(l-) + y/"(l-) = 0. 

The domain T>{A) C C 2 ([0, 1]) is uniquely determined by the boundary con- 
ditions ( |1.1| >, \\.2\ . 

Furthermore the following statement holds true: 



Theorem 1.3. Every Browian motion on I in the sense of definition 1.1 is a Feller 
process. 



Proofs of these theorems can be found in [21 ]. For the case of metric graphs proofs 



of similar statements are provided in |24|. 



We remark that the boundary conditions determined by a Brownian motion in the 
sense of definition are local, because each of the conditions §\A) , ( |1.2| ) only in- 
volves the values of the function / and its derivatives at one or the other endpoints 
of the interval. In fact, this is a consequence of our requirement that the paths of the 
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Brownian motion be continuous up to their lifetime. More general boundary condi- 
tions arise if one allows for jumps from the endpoints back into the interval. They 
have been discussed in the above quoted work by Feller and by Ito-McKean, see 
also (7J|2TJ. 

In order to explain the ideas of how to construct such Brownian motions pathwise 
from a standard Brownian motion B, for the choice I = M + let us first discuss the 
"pure" cases, where only one of the terms in equation ( |1.1| > is present. 

First consider the choice ao = 1, bo = Cq = 0, i.e., the Dirichlet boundary condi- 
tion /(0) = 0. (This condition is actually not permitted according to the statement of 
the theorem, but temporarily we shall consider it nevertheless, see also below.) Then 
this boundary condition can be implemented by killing the standard Brownian motion 
when it reaches the origin, since by our convention /(A) = 0. The reason that this 
case has to be excluded in theorem [L2"l is that due our definition, a Brownian motion 
X on M + is such that X is equivalent to B, and therefore has to be able to reach 0, 
which is not possible for a right continuous process which is killed upon reaching 0. 
Moreover, a process which is killed upon reaching 0, obviously cannot be a Feller 
process. 

The choice ao = 0, 6o = 1, Co = is the Neumann boundary condition /'(0+) = 
0, and it is well-known (e.g., 1 19 p. 40], cf. also section[2]) that this can be obtained 
from a Brownian motion with reflection at 0. 

Finally, it we choose ao = bo = 0, cq = 1, then we find the Wentzell bound- 
ary condition /"(0+) = 0. Obviously, this means that the generator annihilates the 
function at 0, i.e., the semigroup acts trivially there. In terms of the underlying sto- 
chastic process this means that it stops to move there, and so one can implement this 
boundary condition with the Brownian motion B with absorption in 0. 

For the general case (|TTT|) one has to construct a stochastic process which combines 
all these effects. According to (18] Section 2], it was Feller who suggested how to do 
this (at least for the elastic case ao ^ 0, bo / 0, Co = 0): For a reflecting Brownian 
motion one builds in the effects of killing and absorption on the scale of the local 
time at 0. More precisely, for the Brownian motion on M + with all three coefficients 
in ( |1.1| ) non-zero, one uses the local time of the reflecting Brownian motion to define a 
new stochastic time scale which slows down the process when it is at the origin. This 
makes the origin "sticky" for the process, and replaces the full absorption in the "pure" 
case Co = 1. Then, on this new stochastic time scale the process is killed exponentially 
to produce a term involving /(0) like in the "pure" Dirichlet case. These ideas have 
been carried out in (T8j[l9j, cf. also (7][2l). The case of an elastic Brownian motion 
has also been discussed in [20,39]]. 

Before we end this introduction with an overview of the organization of this article 
let us quickly settle one trivial case of a Brownian motion on M + , namely when bo = 
0. The subcase ao = has already been treated above: the Brownian motion B with 
absorption realizes this boundary condition. So consider ao ^ 0, Co ^ 0, and set 
P = ao/co- Define a stochastic process X as follows: Start with a Brownian motion 
at x > 0, and when hitting the origin keep the process there for an independent 
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exponential holding time of rate (3. After expiration of this holding time let the process 
jump to the cemetery point A. It is straightforward to check that this process is 
a Brownian motion on M + in the sense of definition |1.1[ that it has 1/2 times the 
Laplacean as the generator, and that its semigroup U acts at the origin as Utf(0) = 
exp(— /3i)/(0), t > 0. Thus this process implements the desired boundary condition 
0/(0) + l/2/"(0+) = 0. 

The organization of the article is as follows. In sections [2] - [4] we consider the 
interval I = M.+. Elastic Brownian motion (do / 0, bo ^ 0, Co = 0) on M + is 
constructed and analyzed in section [2] In section [3] the Brownian motion with sticky 
origin (ao = 0, &o / 0, cq ^ 0) is considered, and in section|4] we treat the gen- 
eral case where all three coefficients are non-zero. In section [5J following an idea 
of Knight (2TJ, we use the stochastic processes constructed before on the intervals 
(— oo, 1] and [0, +oo) to piece together a stochastic process on [0, 1], which generates 
the boundary conditions (jTTTJ, ( |1.2| >. The most difficult point is the verification of the 
(simple) Markov property of this process, and some parts of the rather lengthy and 
technical proof are deferred to an appendix. There are also appendices about killing 
a process with a perfect homogeneous functional (included here, because we use the 
same method of killing as in 1 20 21 1, for which we need some results which are not 
easily available in the standard literature), about some results related to the Brownian 
local time, and about general heat kernels and their Laplace transforms. 
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2. Elastic Brownian Motion 

Throughout this article we suppose that (f2, A) is a measurable space endowed with 
a family of probability measures P = (P x , x G M), a right continuous filtration F = 
{Ft, t > 0) which is complete with respect to P, and a standard one-dimensional 
Brownian motion B = (B t , t > 0) relative to F, so that P x (Bq = x) = 1 for every 
x G K. Without loss of generality we may assume that (O, A) is equipped with a 
family 9 = (9 t , t > 0) of shift operators for B: B s o 6 t = B s+t for all s, t > 0. In 
particular, our assumptions mean that B is a normal strong Markov process relative 
to F. As usual, we set Foo = cr(Ft, t > 0). 

The local time of B at the origin will be denoted by L = (L t , t > 0), and we 



choose its normalization so that Tanaka's formula (e.g., p0[ p. 205], (34} p. 207], (32 
p. 68]) holds in the form 

d\B t \ = sign(Bt) dB t + dL t , t > 0. 

It will be convenient to assume that both, B and L, have exclusively continuous paths. 
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FIGURE 1 . Simulation of a path of a Brownian motion (blue) and its 
local time at the origin (red). 



Figure [T] shows a simulated path of a Brownian motion on the real line (blue) and 
its local time at the origin (red)|^] All other paths in the figures below will be con- 
structed from this path. Of course, the figures of these paths can only be considered 
as caricatures of the "true" paths. 

We shall often have occasion to use the right continuous pseudo-inverse of L de- 
noted by K: 

(2.1) K r = inf {t > 0, L t > r}, r e R+, 

with the usual convention inf = +oo. The law of K r , r > 0, under Pq is computed 



in appendix |B| (cf . lemma B.l I. Since for fixed r > 0, K r is the moment when the 
local time L increases above level r, and L only grows when B is at the origin, we 
find that Bk t = 0, r G E + . It is not hard to show that the continuity of L entails 

Lk,. = r,r G K+, and for all r, t > 0, 

(2.2) {K r >t} = {L t < r}. 

In particular, the last equation shows that for every r > 0, {K r < t} belongs to 
Tt, and since F is right continuous, K r is an ^-stopping time for all r > 0. For 



All simulations of this article have been done with SciLAB 5.2.1. 
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later purposes we also remark here that due to its right continuity K is a measurable 
stochastic process. 

Let us consider the reflecting Brownian motion \B t \,t > (cf. figure [2]). Clearly, 




FIGURE 2. A path of a Brownian motion (grey) and the associated 
reflected path (blue) 

the transition semigroup of the reflecting Brownian motion has the integral kernel 

(2.3) p N (t,x,y) = p(t,x,y) + p{t,x,-y), t>0,x,y€R + , 
where p(t, x, y) denotes the usual heat kernel on the line, i.e., 

(2.4) p(t, x, y) = —L= e-^-f) 2 / 2 *, t > 0, x, y G K+. 

v27rf 

Hence the resolvent R N = (R^ , A > 0) of the reflecting Brownian motion has the 
integral kernel 

(2.5) rf (x, y) = ^= (V^AI^I + e -V2X(x+y)\ ^ X > 0, x, y £ R + . 

It is now a straightforward calculation to show that for every / G Co(M+) we have 

(2.6) «/)'(0+) = ; A>0, 
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i.e., the generator of reflecting Brownian motion is 1/2 times the second derivative 
with Neumann boundary condition at the origin. 

Elastic Brownian motion B e is the stochastic process defined by killing the reflect- 



ing Brownian exponentially on the scale of the local time at the origin, cf., e.g., 1 19 
p. 45], [2"T] p. 158], pO) p. 425 f.]. Here we shall follow the construction given 
in [20} [21 1, which is slightly different from the one in [19], or in — in the more 



general form of killing on the scale of a perfect continuous additive homogeneous 
functional — [3 , 39j. For a more detailed account the interested reader is also referred 
to appendix [A 

Let /3 > 0, and introduce the auxiliary probability space jB(R+), Pp) where 
P/3 is the exponential law of rate f3. Denote by S the random variable S(s) = 
s, s G M+. Now form for every x G M the product space of (Q,A,P X ) and 
(K+, £>(M+), Pp), and denote the resulting probability space by (O, A, P x ). S and 
all random variables on (Q,A,V X ), x G M, are extended in the trivial way to the 
product space, and we denote the extended random variables by the same symbols. 

On (Cl, A) define a random time by 



(2.7) 

Observe that we can write 



C/3 = inf{t>0, L t >S}. 



(2-8) C/J = K s , 

because K is a measurable process. Now set 



(2.9) 



t < (p, 
t>Cp, 



where A is a cemetery state. In figure [3] the local time of the path of a reflected 



A 




0.8 £ 1.0 



Figure 3 . Construction of elastic Brownian motion 
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Brownian motion (in grey) is drawn in red, the value of S is 0.58 (horizontal dashed 
line). The time of killing the reflected Brownian motion for the corresponding path of 
elastic Brownian motion is at £« = 0.89 (vertical dashed line). The path of the elastic 
Brownian motion is the blue piece of the depicted path, extended to be equal to A 
after time (p. 

As is shown in appendix |a| (Cl, A) can be equipped with a right continuous filtra- 
tion P = (P t , t > 0) (there denoted by (M. t , t > 0)) which is complete with respect 
to the family (P x , x G such that B e is ^"-adapted and a strong Markov pro- 

cess relative to P. In particular, B e is a strong Markov process relative to its natural 
filtration. Moreover, the lifetime (p is a stopping time for P. 

In the sequel we make the convention that every real valued function on M + is 
extended to R + U {A} via /(A) = 0. 

Now we compute the boundary condition of the generator of the elastic Brownian 
motion. Our calculation is based on Dynkin's formula (cf., e.g., |6j p. 133], fl9] 
p. 99], (391 p. 131]) for the generator A e of the elastic Brownian motion B e . 



Theorem 2.1. The domain D(A e ) of the generator A e of the elastic Brownian motion 
B e of parameter j3 > is equal to the space of functions f £ Cq so that 

(2.10) /'(0+) = /3/(0) 
holds. 

Proof. By construction of B e the origin is not an absorbing point for B e . Thus it is 
permissible to apply Dynkin's formula at the origin: for every / G Cq 

(2.11) A^0) = ^ E ^l^- m , 

where r e , e > 0, is the hitting time of the complement of the interval [0, e), i.e., of the 
set [e, +oo) U {A}, by B e . Thus r e = H e A (p, where H e is the hitting time of the 
point e on the real axis by the reflecting Brownian motion \B\. Write 

E (f(B*j) = E {f(B*);H € < + E (f(B*J; H e > Cp) 
= /(e) P (H e < 0?) + /(A) P (H e > Cp) 



where we used lemma B.4 and the convention /(A) = 0. Thus we obtain 



A ' 1 (0) " >SS iik) (lT^ (/(e) " /(0)) " TT* /(0) ) 

40 E {T t )l + ep \ e HJy 1 

We know that the limit on the right hand side exists, and clearly 

_L / /(e) -/(0) _ ^ N = //(Q+) _ ^ 
e4,o i + ep \ e / 
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On the other hand we have -Eo( r e) < Eq(H e ), and the latter expectation is equal to e 2 
(e.g., [ 19, Sect. 1.7, Problem 6], |7) Chap. II, Problems 5, 18], (39} Example III.24.3]), 
so that 

^0(r e ) " e 

Therefore we must have /'(0+) - /3/(0) = 0. □ 

Remark 2.2. In order to compare the boundary condition ( 2.10 ) with ( |1.1| > with Co = 
0, for given do, &o / with ao + bo = 1, we just have to make the choice of (5 above 

as (3 = a /(l - ao). 

Next we compute the resolvent i? e = (i2|, A > 0) and the semigroup U e = 
(Ujr, t > 0) of the elastic Brownian motion. First we prepare with the following 
lemma: 

Lemma 2.3. For all f € Co(K+), A > 0, x £ R+, the following formula holds true 

(2.12) ^/(x) = - e A (x) 4(e- A ^) </(0), 
vv/zere 

(2.13) ex ( x )= Ex {e- XHo )=e-^ x , 

and Hq is the hitting time of the origin by the Brownian motion B. 
Proof. By construction of B e , 

Rlf{x) = E x ( K j\- xt f{\B t \)dt) 

= R%f(x)-E x (J°° e- xt f(\B t \)dt 

POO 

= R^f{x)-E x [e-^ J e~ xt f{\B t+Cp \)dt 
We compute the last expectation 

E x (e- X <e J" ' e- xt f(\B t+Qj \)dt) 

= (3 r e-^E x (e- XK 



f°° \ 
J e- xt f(\B t+Ks \)dt)ds 

J°° e~ xt E x [f[\B t+Ks |) | T K ^j dt) ds 

- XK °)E (J™e- xt f(\B t \)dt) ds 

= P / e-P s E x (e- XK °)dsR%f(0), 
Jo 

where we used relation (2^8 1, the strong Markov property of B with respect to the 
^-stopping time K s , and Bk s = 0. Since L only grows when B is at the origin, we 



r-oo 

(3 J e- ps E x [e 
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have for all s > 0, K s > Hq, whence K s = Hq + K s o 9 Hq . Therefore we get with 
the strong Markov property of B with respect to Hq 



E x (e- XK °) = E X { 



E x {e- XH °)E { 
^ x E (e 



e 



\K 3 



where we used the well-known Laplace transform of the density of Hq under P x , 
e.g., [ 19, p. 26], (20j p. 96] or (34j p. 67]. Inserting the last expression above, we get 
formula pH] ). □ 

Remark 2.4. We want to emphasize that the same calculation works for any Brown- 
ian motion on R + in the sense of definition 1 . 1 with infinite lifetime having a local 
time at the origin, i.e., a PCHAF (cf. appendix \A\ whose Lebesgue-Stieltjes measure 
is carried by the origin. In section [4] we shall take advantage of this fact when with 
the corresponding analogue of formula ( |2.12| ) we compute the resolvent of a Brow- 
nian motion with the boundary condition (jTTTJ in its most general form, i.e., with all 
coefficients non-vanishing. 

To get an explicit expression for the resolvent R\ , it remains to calculate the expec- 



tation Eo(exp(— AC/?))- But this is easy with the help formula ( B.3 1 in appendixiB 



-/3s 



E, 



o(e 



-A A' 



) ds 



2Xs 



ds 



(2.14) 



2A 



Thus we have proved the following 

Corollary 2.5. For all f € C (M+), A > 0, x £ 



(2.15) 



Rif(x) = R^f(x) 



e A (x)</(0) 



P + V2A 



With formula ( |2.15| > it is now very easy to give another proof of theorem |2Tj Let 
G C (M+), A > 0. Then R% f satis 
therefore we get with equation ( |2.13 1 



/ G Co(K+), A > 0. Then R^f satisfies the Neumann boundary condition at 0, and 



(2.16) {R e xf)'(0+) 
while on the other hand 

(2.17) R e x f(0) = 
and therefore 



/3V2A 



</(0), 



P + V2A 



2A 



«S7(o), 



(2.18) 



P + V2A 
R e x f)\0+) = PR e x f(0). 
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Since for every A > 0, R e x maps Cb(M+) onto the domain of the generator of B e , we 
have once again proved theorem |2~Tj 

From equation ( |2.15| > we can read off the resolvent kernel rf (x, y), x, y G M.+, of 
the elastic Brownian motion 

J? (a- ,.\ _ 2/3 _L p -V2X(x+y) 



(2.19) r e x (x,y) = r%(x,y) - " -= e -v^*+fl x, y G R+, 

p + v 2A \ 2A 

where r¥(x, y) is the Neumann kernel (2"3l. In order to compute the transition kernel 
of the elastic Brownian motion, we want to find the inverse Laplace transform of 
r\{x, y) as a function of A > 0. It turns out that it is more convenient to achieve this 
if we rewrite rf(x, y) in terms of the Dirichlet kernel 

(2.20) rf (x,y) = -L ( e -^l*-i/l _ e-v^A^+w)^ x, y £ R + . 

v 2A 

Then for x, y G R+, 

(2.21) r£(s >y ) = rf (x,y) + - 2 e -^+y), 

p + v 2A 

Now the inverse Laplace transform of the second term on the right hand side is given 
by formula (5.6.12) in p). So we find for the transition density p e (t, x, y), t > 0, x, 
y G M+, of the elastic Brownian motion 

(2.22) p e (t, x, y) = p D (t, x, y) + 2g pfl (t, x + y), 

where p (t, x, y) is the transition kernel of a Brownian motion on M + which is killed 
when reaching the origin, i.e., whose generator is one half times the Laplacean on M + 
with Dirichlet boundary conditions at the origin: 

(2.23) p D (t, x, y) = p(t, x, y) - p(t, x, -y). 
Moreover, we have set 

j3 / P 2 t\ fx ft\ 

(2.24) gp,o(t,x) = g{t,x) - - exp^x + —J erfc^-^= + fiyj - J, 

and g(t, x) denotes the standard GauB-kernel 

(2.25) g(t, x) = —L= e~ x2/2t , t > 0, x G M. 
Alternatively, we can write p e (t, x, y) in the form 

(2.26) p*(t, x, y) = p N (t, x, y) - erfc + (jjt) , 



where p (t, x, y) is the Neumann heat kernel (2.3 1. 
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3. Brownian Motion With Sticky Boundary 

In this section we consider the case ao = 0, bo, Co > in the boundary condi- 
tion As mentioned in the introduction, this case is handled by introducing a 
new time scale which slows down the reflecting Brownian motion \B\ when it is at 
the origin. Again, the idea of how to implement this with the help of Levy's local time 
goes back to the fundamental paper [ 18] of Ito and McKean. We continue to denote 
the local time of the Brownian motion B at the origin by L = (L tl t > 0). Let 7 > 0, 
and set 

(3.1) r-\t) = t + 1 L u t>0. 

Since L is increasing, r _1 is strictly increasing. Furthermore, t _1 (0) = and 
t _1 (+oo) = +00 are valid, which implies that r exists and is strictly increasing, 
too. In figure [4]the local time of the path of a reflected Brownian motion from above 
is depicted in red, the path of t _1 (with 7 = 0.3) in green, while the path of the new 
time scale r is in blue. (The grey diagonal line represents the "deterministic time 
scale".) 




FIGURE 4. Construction of the time scale r (blue) from the local 
time L (red) and r _1 (grey) 

Define a Brownian motion B s on R + with sticky boundary at the origin by 
(3.2) B s t =\B T{t) l t>0. 

Figure [5] shows the path of a reflected Brownian motion from above in grey, and 
the corresponding path of a Brownian motion with sticky boundary at in blue. By 
formula ( |3.1[ ) it is clear that r _1 increases at a faster rate than the "deterministic time 
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0.5 




FIGURE 5. A path of a "sticky" Brownian motion (blue) constructed 
from the path of a reflecting Brownian motion (grey). The local time 
at of the "sticky" path is in red. 



scale t h-» t" whenever B is at the origin, and therefore r increases slower than the 
deterministic time scale in these instances. Hence — at least in a heuristic sense — 
B s experiences a slow down when it is at the origin. Below this remark will be made 
precise. With these considerations in mind, we shall also call 7 the parameter of 
stickiness of B s . 

First we investigate the (strong) Markov property of B s , and we fill in some details 
of a proof given in f2~T] p. 160]. We recall the additivity property 

(3.3) L s+t = L s + Lt o 9 S , s,t£R + , 



of the local time L, which directly leads to the following formula, see [21 p. 160], 



r(t) 



7 r(s)> 



s, t > 0, 



(3.4) r(s + t) = r(s) 

i.e., r is additive on its own scale. 

For t > define F^ = T T uy Then we have the following result 

Lemma 3.1. For every t > 0, r(t) is an J 7 — stopping time. Furthermore, F s = 
(J~i , t > 0) is a filtration of sub— a— algebras of J- 00, which is right continuous and 
complete relative to P, and such that for every t > 0, J 7 / C J~t- Furthermore, 

Proof. First we show that for every t > 0, r(i) is an ^-stopping time: let s > 0, then 
by dtr) 

Mi) <s} = {t< t- x (s)} = {t < s + jL s }, 

and the last set belongs to F s . In particular, P T (t) is a sub-c-algebra of F^, and 
this entails the last statement of the lemma. That F s is a filtration follows from the 
fact that t is increasing, and that for every t > 0, J 7 ^ C is due to r(t) < t. 
Moreover, r(0) = gives Fq = Fq, and therefore the completeness of F entails 
the completeness of F s . It remains to show that F s is right continuous. To this end 
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let t > 0, and assume that A G Ff + , i.e., A G ^neN^ 7 ^!^- We have to prove that 
A G J 7 / = .?>(*)• The fact that r is (pathwise) continuous implies that r(t + 1/n) 
converges from above to r(i) as n tends to infinity. Hence for u > we get 

{r(t)<u}= ]J{T(t+l/n) <u}= |J J{r(t+l/n) < u - e}, 

neN e>0, eSQnSN 

and 

An{r(t)<n}= [J y An{r(t + l/n) <u-e}. 

<=>0, eeQnSN 

Because A G Jv(t+i/n) for a11 ^ G N, we find that for all n G N, e G Q with e > 0, 

A n {r(t + 1/n) < u - e} G F u - t C 

Thus we have shown that for all u > 0, A n {r{t) < u} G F u . But then for all u > 0, 
A n {r(t) < u} G J r «+. By assumption = J 7 ^, and this gives A G TrU)- ^ 

Since B has continuous paths, B is progressively measurable with respect to J 7 , 
and therefore (see, e.g., |34| Proposition 1.4.9]), for every t > 0, |-B T ( t ) | is measurable 
with respect to F r ^, i.e., B s is adapted to F s . Henceforth we shall always consider 
B s relative to the filtration F s , unless mentioned otherwise. 

Next we define a family 6 s = (8f, t > 0) of mappings from into itself by 

(3.5) s t=Or(t), t>0. 

Then we obtain from equation (^l for s,t£ R + , 

B t+t = \ B r( S )+r(t)oe T{s) \ 

(3-6) = |B T(t )| o T ( s) 

and therefore 

(3.7) B s s+t = B s t o9l 

which shows that 6 s is a family of shift operators for B s . 

Now we claim that B s is a Markov process relative to F s . First we remark that it is 
an easy exercise to show that the reflecting Brownian motion |J3| is a strong Markov 



process with respect to T . Let x G M+, s, t > 0, C G B(K+). Recall lemma 3.1 



by which r(t) is an ^-stopping time for every t > 0. Using equation (3^)1 we can 
calculate as follows: 

P x {B s s+t G C | T t s ) = P X {\B T(S) \ o 9 T(t) G C | F T(t) ) 

= P\B T{t) \{\B T{s) \eC) 

= P BI (B* a €C), 

and our claim is proved. 

Recall that Hq is the hitting time of the origin by the standard Brownian motion B. 
As a next step we want to prove that B s has a strong Markov property with respect to 
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the hitting time Hq of the origin and the filtration T, i.e., we want to show that for all 

x G R+, t > 0, C G B(R+), 

(3.8) P x (B s t+Ho G C | T Ho ) = Po {Bt G C) , P x -a.s. on {H < +oo}. 

All considerations below are done on the set {Hq < +00} which has i^-measure 
one for all x G M+. Observe that t{Hq) = r~ 1 (f^o) = Hq which follows from 

t- 1 (H ) = H + 1 L Ho =H q , 

because Ljj = 0. In particular, it follows that Hq is the hitting time of the origin for 
both B and B s . Thus by (53) 

= l 5 r(// )+r(t)oe T(H()) l 

= l 5 // +r(t)oe Ho l 

= |^T(t)| & H 

(3.9) =B s t o9 Ho - 

Using the fact shown above that for t > 0, S| is J^-measurable, and the strong 
Markov property of the reflecting Brownian motion |J3|, we now compute for x G M+, 
C G B(K+), t > 0, as follows: 

Px [B s t+Ho EC\T Ho )=P x {Bf o Ho G C I Fh ) 

= Px{\B T (t)\ °H I ^ff ) 
= P (|S r(t) |GC) 

= P ( J B t s eC). 

Since Hq is also the hitting time of by B s , it is also a stopping time for F s . Thus 
Ffj Q is a well-defined cr-algebra, and A G J-j^ if and only if A G J 7 ^, and for all 
t > 0, A n {i7 < G J 7 /. Since J 7 ^ C Too (s.a.), we immediately get A G J 7 ^,, 
i.e., J 7 ^ C Fh - Consequently B s is a strong Markov process relative to Hq and the 
filtration F s , too. 

Using the facts established above that {B S ,F S ) is Markovian, and strongly Mar- 
kovian for the hitting time of the origin, we are now in position to apply the same 
arguments as in the proof of Theorem 6.1 in (2T| (cf. also the proof of theorem 4.3 
in (24) with a slightly different argument for the case of metric graphs and somewhat 
more details than those given in [21]) to conclude that B s is a Feller process. But this 
proves 

Theorem 3.2. B s is a normal strong Markov process with respect to F s . 
The necessary path properties being obvious, we have also proved 



Corollary 3.3. B s is a Brownian motion on M + in the sense of Definition 



1.1 
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Next we turn to calculate the boundary conditions of the generator A s of B s . Again 
we will use Dynkin's formula. We first argue that the origin is not a trap for B s : Let 
Z C M+ denote the set of zeroes of the underlying Brownian motion B. Given s > 0, 
choose to > s in the complement Z c of Z. Consider t = r _1 (to)» i-e., t = to + jL to . 
Obviously t > s, and r(t) G Z c . Therefore Bf = \B T u\\ / 0. Hence is indeed not 
a trap for B s . For e > let H^ j6 denote the hitting time of the complement of [0, e) 
by B s , that is, since by construction B s has infinite lifetime, H l t is the hitting time 
of the point e G M + by B s . Then for / 6 Cn(M+) we can calculate as follows: 



^ ^ 40 Eq{H 1 ^ 



lim 



^o(/(^j )-/(0) 

/(e) - 7(0) 



40 Eo(Hy !e ) 
Now we use the following 

Lemma 3.4. Let H e denote the hitting time of e > by the reflecting Brownian 
motion \B\. Then under Pq, 

(3.10) H li€ = H e + 7 L He . 

Proof. By definition 

H j,e = inf{s > 0, \B T ^\ = e}, 
H e = in£{t > 0, \B t \ = e}. 

Set s = H e + jLh c , then r(s) = i.e., |-B T ( S )| = e, and therefore we get the 
inequality 

ff 7) e < H e + jL He . 
Now assume that the last inequality is strict. This would entail 

r(ff 7 , £ )<r(ff £ + 7 L fft )=ff e 

because r is strictly increasing. But \B T r H e )| = e, in contradiction to the definition 
of H e as the hitting time of e by \B\. Thus formula ( 3. 10| > is proved. □ 



We already know (cf. the proof of theorem 2.1 ) that Eq(H e ) = e . On the other 



hand, lemma B.2 states that under Pq, Lh^ is exponentially distributed with mean e. 
Therefore we find by lemma |3~4[ 

(3.11) Eo(H^ t ) =e 2 + 7 e. 

Thus 

(3-12) 2 / "( 0+ ) = 7 ^/(0) = /'(0+), 

and we have proved 

Theorem 3.5. The domain T>{A S ) of the generator A s of the Brownian motion B s 
with sticky origin and stickiness parameter 7 > is equal to the space of functions 
f £ Cg(M+) so that the boundary condition (3.12 1 holds true. 



18 



V. KOSTRYKIN, J. POTTHOFF, AND R. SCHRADER 



Remark 3.6. Assume that in §1.1) we are given ao = 0, &o> Co 7^ 0. Then we choose 
7 = Co/ (1 — Co), and the so constructed process B s has a generator with a boundary 
condition given by \\A) . 

Remark 3.7. Equation ( |3.1 1| > shows that when starting at the origin the Brownian 
motion B s in the average needs the time e 2 + 76 to leave the interval [0, e), while the 
reflecting Brownian motion needs the time e 2 . On the other hand, consider the case 
that B s starts at x > 0, and let e > be small enough so that e < x. Then on the 
average B s takes the time e 2 to leave the interval (x — e, x + e), because there it is 
equivalent to \B\. This shows very clearly the effect of "stickiness" of the origin for 
B s : Not only takes B s longer than to leave a neighborhood of the origin, even it 
does this on a time scale different from that of \B\. 

Finally we compute the resolvent R s = (R x , A > 0) and the transition kernel 
p s (t, x, y), t > 0, x, y G M+, of B s . It will be convenient to introduce for measurable 
functions /, g on M + the inner product 

roc 

(9 J) = / 9(x)f(x)dx, 
Jo 

whenever the integral on the right hand side is defined. 
Lemma 3.8. For all f G Cb(M+), the relation 

(3.13) \ (# A /)"(0+) = -=± (2A (e x , /) - V^A/(0)), A > 0, 

* v 2A + 7A 

holds true, where e\ is given by equation ( |2.13| ). 

Proof. The proof given here is basically the proof in ('21], though somewhat reorga- 
nized and simplified. Let A s denote the generator of B s . Since for A > 0, R s x f 6 
V(A S ), A S R S J belongs to C (M+). Also for x > 0, A s R s x f(x) = 1/2 (R{f)"(x). 
Thus 1/2 (R 8 x f)"(0+) = A s R s x f(0), and with A S R S X = XR S X - id we can compute 
as follows 

l(^/)"(0+) = A J Ri/(0)-/(0) 

= \E (J\- Xt (f(B s t )-f(0))dt) 

= XE (^°° e -*(.+7£.) (/(|S S |) - /(0)) (cfa + -ydLs)) 

= A^o(^°°e- A (^)(/(| J B s |) - /(0)) ds). 

In the last step we made use of the fact that L only grows when \B\ is at the origin. 
Therefore the terms with f(\B s \) and /(0) cancel each other in the integral with 
respect to dL s . Now we use the well-known joint law of (\B S \, L s ) under Pq 

(3.14) P (\B S \ £ dx,L s edy) = 2^H=e-^^ 2s dxdy, x,y>0, 

V2irs 6 
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which follows directly from Levy's theorem (see, e.g., pO) Theorem 3.6.17]), and 
the joint law of a standard Brownian motion and its running maximum under Po 
(e.g., [20, Equation 2.8.2]). Thus we get 

-i roc roo roo 

- (R{f)"(0+) = 2A / / / e-^) (/(*) - /(0)) 
z Jo Jo Jo 

We make use of the well-known Laplace transform (e.g., (33} 1.5.30]) 

(3.15) / e - Xs ^^e~ a2 / 2s ds = e~ V ^ a , a > 0, 
and get 

i roc 

- (Rif)"(0+) = 2A J e~ A ™ e-^(^) (/(*) - /(0)) 

= - (2A(e A ,/)-v / 2A/(0)), 

V2A + 7A 

which proves the lemma. □ 
Corollary 3.9. For all f £ C (M+), A > 0, 

(3.16) R{f(0) = 1 (2(e A ,/) +7/(0)). 

V 2A + 7A 

Proof. Formula ( |3.16 l follows directly from equation ( 3.13 1, the identity Rif = 
A _1 (/ + A s R s x f), and the fact that A s acts as one half the second derivative on 
its domain. □ 

Now we apply the well-known first passage time formula (e.g., fl9"l p. 94]) in the 
following form (for the convenience of the reader we include a detailed proof of this 
formula within a general context in appendix |D]| 

(3.17) R{f(x) = fl?/(s) + E x (e- XHo ) R{f(0), 
which combined with ( |2.13[ ) simplifies to 

(3.18) R{f(x) = fl?/(s) + e x (x) R s J(0). 

From this formula we can read off the resolvent kernel r s x (x,dy) of the Brownian 
motion with sticky origin: 

Corollary 3.10. For all A > 0, x, y € R+, 

(3.19) r{(x,dy) = r°(x,y)dy+ 1 - (2e~^^ + je~^ e (dy)), 

V 2A + 7A 



where rP(x, y) is defined in (2.20), and eq is the Dirac measure in 0. 
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From the first passage time formula ( |3.18| >, we immediately get for all / G Cq( 



(3.20) {R s xf)'(0+) = 2(e A , /) - V2XR a J(0), A > 0. 

Inserting formula ( |3.16[ ) into the right hand side and comparing the result with ( 3.13 1, 
we find 

(3-21) («i/)'(0+) = |( J Ri/) ,/ (0+), 



which is another proof of theorem [33] 

The inverse Laplace transform of the second and the third term on the right hand 
side of equation ( |3.19| ) can be done with the help of formula (5.6. 16) in [8]. The result 
is the following 

Corollary 3.11. For t > 0, x, y E M+, the transition kernel p s (t,x,dy) of the 
Brownian motion with sticky origin is given by 

(3.22) p s (t, x, dy) = p D (t, x, y) dy + 250,7 (t, x + y)dy + 750,7 %) eo(dy), 
where p D (t, x, y) is given in ( |2.23[ ), and 

„ . , 1 /2x 2t\ r / x y/2t \ 

(3.23) go~(t,x) = - exp 1 — 5- eric — = H , t > 0, x > 0. 

7 v 7 7 ' v V 2t 7 ' 

Remark 3.12. #0,7 is a special case of a more general heat kernel which will play an 
essential role in the next section, and which is investigated in detail in appendix [C] In 
particular, it is shown there, that as 7 J, 0, 50,7(^7 x ) converges pointwise to the usual 
GauB kernel g(t, x) ( |2.25[ ). Therefore, as 7 \. 0, the transition kernel of the Brownian 



motion with sticky origin converges pointwise to the Neumann transition kernel, i.e., 
to the transition kernel of the reflecting Brownian motion, as it should. 

We conclude this section with a discussion of the local time of {B S ,F S ) at the 
origin, and this will also provide results needed for the construction of the general 
Brownian motion on M + in the next section. Define L s = (Lf, t > 0) with 

(3.24) L s t =L T{t) . 

Observe that L s has continuous non-decreasing paths and that Lq = 0, because 
t(0) = and Lq = 0. Moreover, L s is ^-adapted since for every t > 0, r(t) 
is an ^-stopping time, so that the progressive measurability of L relative to T entails 
that L T u\ is measurable with respect to F T (t) = ^i- W e claim that L s is additive with 
respect to the shift operators 9 s . Namely, L s s+t = L s s + Lf o Q s s holds for all s, t £ M + , 
as a consequence of the relations ( |3.3| >, ( |3.4| ), and the definition 9% = 9 T uy Indeed, 
by evaluating at uj G Q we obtain 

K+tiu) = L T ( s+ t)(a;)(a;) 

= ir(s)H+r(t)(e T(s)M M)(w) 

= I- T ( a )( w )(w) + £ T(t)(0 T( , )(w) (w))(^TW(a;)(w)) 

= LfH + Lf(^H). 
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Altogether we have proved that L s is a PCHAF for (B s , J 73 ). Finally we show that 
the Lebesgue-Stieltjes measures of its paths are carried by the origin, that is, we have 
to prove that L s only increases when B s is at the origin. But L s increases at time 
t > if and only if L increases at time r(t) which is true only if |-B T (j)| = Bf is at 
the origin. Hence L s is a local time at for B s . We determine its normalization by 
computing its a-potential (see e.g., [3]): Let a > 0, then 



Er 



-at 



dL s t 



Eq 
Eq 



e- at dL T{t) 



-at—a"fLt 







1 



= E 

07 

= Eq 

aj 7 

The well-known law of L t under Pq is given by 



(3.25) 



Po(L t £ dy) 



'2id 



-y 2 /2t 



dy, 



dh tj 

e - at -a-,L t d , j 

y >o, 



as, for example, we can easily deduce from the joint law ( |3.14| ). Thus we can compute 
the above expectation 



E ( 



-at—arfLt 



dt 



-at 



g(t,y)dt) dy 



V2a (v2a + arj) 
where g{t, y) is the GauB kernel ( 2.25 1. Hence 

1 



-at 



dL) 



x G R + , a > 0. 



'2a + cry 

Now consider x > 0. Recall that above we had argued that on [0, Hq] the paths of B 
and B s coincide, so that Hq is also the hitting time of B s , and L s is identically zero 
on [0, Hq], too. Hence the strong Markov property of B s gives 



E, 



-at 



dLl) =EJe 



-aHo 



Eje~ aHo E, 



-at 



dL 



E x (e- aHo )E 



-at 



t+H Q 



dLf) . 



H 



Thus we finally obtain 



(3.26) 



E, 



-at 



dL s t 



2ax 



/ 2a + 07 ' 



x G R + , a > 0. 
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4. The General Brownian Motion on 



In this section we construct and investigate the most general Brownian motion on 
M + in the sense that all three coefficients in the boundary condition ( fLl) are non- 
zero. That is, in contradistinction to our discussion in the previous section now the 
term /(0) appears, and we have already seen in section [2] that such a term appears 
when killing the process at the origin on the scale of the local time at 0. We proceed 
as follows. 

Fix 7 > and consider the Brownian motion B s defined in section [3] Fix f3 > 0, 
and consider the product spaces (Cl, A, (P x , x 6 , and the random variable S as 
in section [2] As there, extend the random variables Bf,Lf,t> 0, to Q in the trivial 
way. Set 



(4.1) 
and 



C A7 = mf{t>0, L\ > S}, 



(4.2) 



if t < C/3,j, 

if* > C/3, 7 - 



Figure [6] shows this construction on the basis of the path of a "sticky" Brownian 
motion depicted in figure[5](which in turn was constructed from the path of a reflecting 
Brownian motion in figure [2]). In figure [6] the local time L s at of the motion B s is 
in red, its path in blue and grey. When L s increases above level S (given in the figure 
by S = 0.58) at time 7 = 1.06, the process is killed, i.e., we obtain the blue piece 
of path, and after time Qp a = 1.06 the path is constant and equal to A . 



A • 




FIGURE 6. A path of the general Brownian motion on IR + (blue) 
constructed from the path of a sticky Brownian motion (grey) and its 
local time (red) 
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Again we can apply the results in appendix [X| In particular, we can equip (Cl, A) 
with a right continuous, complete filtration F 9 which is constructed there, so that 
B 9 is adapted and strongly Markovian relative to F 9 . Moreover, B 9 has obviously 
the right path properties, and it is equivalent to a standard Brownian motion up to 
its hitting time of the origin. Thus B 9 is a Brownian motion on M + in the sense of 
definition ll.il 

Next we compute the boundary conditions of the generator A 9 of B 9 . In this case 
it seems to be difficult to do this via Dynkin's formula. Therefore this time we only 
give an argument based on the computation on the resolvent R? of B 9 . However, 



instead of using the first passage time formula as is done in |T8j[T9 21 1, here we use 
the analogue of formula ( |2.12| ), which in the case at hand reads 

(4.3) R{f{x) = R{f(x) - e x (x) E {e~ x ^) R{f(0), 

where R s is the resolvent of B s . As mentioned in remark |2T4{ the proof is the same 
as for the elastic case discussed in section [2] Thus we only need to compute the 



expectation occuring in equation ( |4.3| ). To this end we first prove the following result. 

Lemma 4.1. For all /3, 7 > 0, 

(4-4) C/3, 7 = C/3 + 7S 

holds true. 

Proof. Define the subsets 

M S ) = {t>0,L t > S} 
Jp,v(S) = {t>o, v t > S} 

of M.+. Since U t = L r (t) and r is a bijection from M + onto itself (cf. section [3j>, we 



get Jp{S) = if and only if Jp tl (S) = 0, and in this case relation (441 is trivially 
valid due to the convention inf = +00. Now assume that one of the sets is non- 
empty, and therefore both are non-empty. Then we have t € J/3^(S) if and only if 
r(t) G Jp(S), and therefore r is also a bijection from J^^(S) onto Jp{S). Because 
r is strictly increasing (cf. section [3]) it follows that 

r(inf J/3,~/(S)) = inf Jp(S), 

or in other words 

The continuity of L implies = S, and the proof is concluded. □ 
Corollary 4.2. For all A > 0, 

(4.5) E (e- X ^) =MA), 
holds true, where p is given by 

(4.6) p(A) = (/3 + v / 2A + 7 A)" 1 , A > 0. 
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Proof. By formula (2.8 1 and lemma |4~T 



poo 

= f3 e ~(?+^> E (e- XK °) ds 
Jo 

= /3 [°° e- [P+V ^ +lX)s ds 
P 



o 



/3 + ^/2A + 7 A' 

where again we used equation ( |B.3| ). □ 

Equation ( |43] > immediately yields for / G Cq(R+), A > 0, x G R+, 

(4.7) R{f{x) = R{f(x) - f3p(X)e x (x)R{f(0), 
and in particular for x = 

(4.8) R{f(0) = (V2A + 7 A)p(A)^/(0). 
Differentiation of the right hand side of equation ( |4.7[ ) at from the right gives 

(4.9) ( J Ri/)'(0+) = [R s x f)'(0+) + f3V2Xp(X)R s J(0). 

In order to compute the second derivative of R g x f at from the right, we make use of 
the fact that R s x f satisfies the boundary condition ( |3.12[ ): 

= {Rlf)'(0+)- 1 \Pp(\)Rlf(0) 

(4.10) = {R{f)'(0+) + f3V2\p(\)R{f(0) - P(V2X + 7 A)p(A)i?|/(0). 
Upon comparison of equation ( |4.10| ) with ( |4.8| ) and ( |4.9[ ), we find the following result. 

Theorem 4.3. The domain D(A 9 ) of the generator A 9 of the Brownian motion B 9 
with parameters f3, 7 > is equal to the space of functions f G Cg (M+) so f/zatf 

(4.11) / g/( )-/'(o+) + |/"(0+) = 0. 

Remark 4.4. In order to compare the boundary condition ( 4.11 1 with ( |1.1| >, for given 



non-zero ao, 6o> cq, as in theorem 1.2 we only need to choose (3 = ao/bo, and 



7 = cq / 60 . to get the boundary condition in the form ( |1.1| >. 



As another byproduct of equation ( |4.7| ) we get in combination with formulae ( |3.16[ ), 
( |3.18| ) the following result. 

Corollary 4.5. For all A > 0, / G C (M+), 

(4.12) iif/(x) = i?f/(x)+p(A)(2(e A ,/)+ 7 /(0))e A (x), x G M+, 
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and in particular, R 9 X has the integral kernel 



r{{x, dy) = rf (x, y) dy + 2p{\) e A (x + y) dy 



(4.13) 



+ lpWe\(x)e (dy), x, y G K+, 
where p(X) is given in (4.6 1, e\ in (2.13), and in ( |2.20[ ). 

In order to determine the transition kernel of B 9 we have to compute the inverse 



Laplace transform of the right hand side of equation ( |4.13| >. This can be done with 
formulae provided and proved in appendix [C] For t > and x G M + set 



(4.14) 9f3,j(t,x) 



i i r s + 'jx 



7 2 Jo (t~ s) 3/2 



exp 



(s + 7x) 2 
'2 7 2 (i-s) 



-? s h ds, 



Then we get from corollary C.2 



Corollary 4.6. For t > 0, x, y G M + , f/ie transition kernel of B 9 is given by 
(4.15) x, dy) = p D (t, x, y) dy + 2# /3>7 (i, x + y) dy + 7^,7 (*> x ) e o(dy)- 

Remark 4.7. In appendix [C] it is argued that for all i > 0, x G M+, 

limffft^M) = 50,7 (*,^) 



lim^ i7 (t,x) = gp,o(t,x) 

74-0 

limlim o« ^(t, x) = limlim qa ~(t, x) 
Pio 740 1 n 740 /340 p ' 7 



g(t,x), 



where gpfl{t, x) is defined in ( |2.24| ), go,~/(t, x) in ( |3.23[ ), and x) in ( |2.25| >. Thus, 
in the limit that we let one or both parameters (3 and 7 converge to zero from above, 
we recover the Brownian motions discussed previously, as we should. 

5. The General Brownian Motion on a Finite Interval 

In this section we construct a Brownian motion on a finite interval following an 
idea of Knight sketched in (2TJ. There will be no loss of generality to choose the 
interval [0, 1]. From the discussion in the previous sections we may consider as given 
the following building blocks 

(a) a standard one-dimensional Brownian motion O = (Ot, t G 

(b) a Feller Brownian motion A = (A t , t G M+) on [0, +00) which implements 
the boundary conditions (jXTTJ) at 0, 

(c) and a Feller Brownian motion B = (B t , t G M+) on (— 00, 1], which imple- 
ments the boundary conditions ( |1.2[ ) at 1. 

Roughly speaking, the idea of construction is as follows. The process is obtained 
as the standard Brownian motion O starting at some point x G [0,1], and it runs 
until it hits one of the boundary points or 1, say 0. From that moment on, the 
process continues as the Feller Brownian motion A. If A hits 1 before its lifetime 
expires, it continues from the moment of hitting 1 as the Feller Brownian motion B. 
If this process hits before its lifetime expires, another independent copy of the Feller 
Brownian motion A is started at this moment in 0, and so on. 
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This construction has similarity with some aspects involved in Ito's excursion the- 
ory p7| , cf. also 1 35]- 37j. But the situation considered here is much simpler, and 



can be handled with relatively elementary — even though somewhat lengthy — argu- 
ments. 

In the next subsection the construction sketched above is formalized. In subsec- 
tion [5]2] it is proved that the resulting process has the simple Markov property, some 



technical details being deferred to appendix \E\ The strong Markov property is proved 



in subsection 5.3 where we also compute the generator of the process. 



5.1. Construction. Let I denote either of the intervals (— oo, 1], [0, +oo) or [0, 1], 
and write Ca(K + , /) for the space of mappings / from M + into / U {A} which are 
right continuous, have left limits in I, are continuous up to their lifetime 

C/ = inf{t>0, f(t) = A}, 

and which are such that f(t) = A implies /(s) = A for all s > t. If / = [0, 1] we 
shall also simply write Ca (!&+)■ 

Assume we are given a family (M,M,p = (p x , x G [0, 1])) of probability spaces, 
on which a standard one-dimensional Brownian motion O = (Ot, t G M + ) with ex- 
clusively continuous paths is defined. Moreover, we assume that A = (A t , t G M.+) 
and B = (B t , t G M+) are two Feller Brownian motions on [0, +oo), (— oo, 1] 
respectively, constructed from O as in the previous sections, and implementing the 
boundary conditions ( |1.1| >, ( |1.2| >, at 0, 1 respectively. It follows from the construc- 
tion that A and B have all their paths in C A (M + , [0, +oo)) and C A (M+, (-oo, 1]) 
respectively. 

The hitting time of c = or 1 by O will be denoted by p c , and we set p = po A p\. 
We write a for the hitting time of 1 by A, and r for the hitting time of by B. 

Below it will be convenient to have somewhat specialized versions of these pro- 
cesses A and B at our disposal, and we construct these now. Let (H^C*, Ql), c = 0, 
1, be two probability spaces. Furthermore, let A 1 be a Feller Brownian motion on 
(Hq,Cq,Qo)> with all paths belonging to CaQR+, [0, +oo)), starting in 0, and such 
that under Q\, A 1 is equivalent to A under ^o- Analogously, B 1 is a Feller Brownian 
motion on (E\, C\, Q\), which under Q\ is equivalent to B under pi, and which has 
exclusively paths starting in 1 and belonging to Ca(K+, (— oo, 1]). Let (H 1 , C 1 , Q 1 ) 
denote the product probability space, and view A 1 , and B 1 as well as their hitting 
times a 1 , r 1 , of 1, respectively, as random variables on the product space. 

Let 

(3°, C°, Q° = (Q° x ,x G [0, 1]), O , pi pi p°) 



be a copy of 
and let 



(M,M,p = (p x , x G [0, l}),0, P o, Pl ,p), 



({E k ,C k ,Q k ,A k ,a k ,B k ,T k ), k G N) 
be a sequence of copies of 
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Define the product space (H, C) as 

H = X C = (g) C k . 

Furthermore, on (H,C) we consider the family Q = (Q x , x G [0, 1]) of product 
probability measures defined by 

(5.1) Q x = Q x ®((g)Q k ). 

ken 

All stochastic processes and random variables introduced above on the measurable 
spaces (E k ,C k ), k G No, are also viewed as defined on (H,C). 
Let us set 

(5-2) H = {p° < pi}, Hi = {p\ < p° }, 

and observe that both sets belong to C. (We remark that in the second event equality 
occurs if and only if a path of O never hits both endpoints of the interval, i.e., if 
Pq = Pi = +oo. Of course, these paths form a negligible event.) 

Next we define on (H, C, Q) a sequence (Sk, k G No) of random variables with 
values in [0, +oo] which will serve as crossover times for the pieces of paths from 
which we shall construct the first version of the stochastic process we aim at. We set 
So = 0, Si = p°, and for k G N, uj G H, 

a k (uj), uj G Ho, k odd, or uj G Hi, k even, 

r k (uj), otherwise. 



(5.3) S k+ i(ui) = S k (uj) 



By an application of the Borel-Cantelli-Lemma, for example, it is easy to see that 
there exists a set H' G C, so that for every x G [0, 1], Q x (z!) = 1, and for all uj G H', 
the sequence (Sk(uj), k G No) is strictly increasing to +oo. 

Now we can construct a "raw model" Y = (Y t , t G M+) of the stochastic process 
desired. On H \ H' we define Y% = A for all t G M+. Given t G M+ and ui G H', there 
exists /c G N so that t G [Sk(w), S k+ i(uj)). Then we set Y t (u) = 0°(w) if k = 0, 
and if k G N, 



(5.4) Y t (cj) 



A k _ Sk (uj), uj G Ho, k odd, or uj g Hi, k even, 
^B k _ Sk (u), otherwise. 



Moreover, we make the convention Y +OQ = A. (Note that due to their path properties, 
A k and B k are measurable stochastic processes, so that A k _ Sk and B k _ Sk are well 
defined measurable mappings, and hence Y t is indeed for every t > a random 
variable on (H, C) with values in [0, 1] U {A}.) Figure [7] schematically explains the 
idea: The piece of the (unrealistic) trajectory in green at the beginning belongs to 
O , the blue pieces to trajectories of copies of A, and the red pieces to trajectories of 
copies of B. 

We want to argue that once Y has reached the cemetery point A, it cannot leave it. 
To this end, let us assume that uj G Ho, that k is odd, and that t G [Sk(uj), S k +i(uj)). 
Then Y t {uj) = A k _ s (uj). Suppose that the lifetime of A k (u) expires before this 
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Figure 7. Construction of Y 



trajectory hits 1, and hence Y s (uj) = A for some s in the interval [Sk(uj), Sk+i(w)). 
Then a k (uj) = +oo, and therefore also Sf.+x(oj) = +oo. Thus there are no more 
finite crossover times for this u, and hence after s the trajectory Y. (uj) is identical to 
the trajectory of A k _ Sk (uj), namely constantly equal to A, as it should. It is clear that 
in all the other cases one can argue in the same way. 

Observe that by construction all paths of Y belong to Ca(R+), and in particular 
they are all right continuous. Hence Y is a measurable stochastic process, and it 
follows that if R is any random time (i.e., a random variable on (S, C) with values in 
M + U {+oo}) the stopped process Y. Aj r = (Y^r, t > 0) is a well defined stochastic 
process. 

Consider the measurable space (fl,A), where U = Ca(K + ), and A is the a- 
algebra generated by the cylinder sets of Ca(R+). Clearly, Y, considered as a map- 
ping from S into Q = Ca(M+) is C/^4-measurable. We denote by P x , x £ [0,1], the 
image of Q x under Y. Furthermore, by X = (Xt, t G M+) we denote the canonical 
coordinate process on (7a(K+). We make the convention X +OQ = A. 

Assume that we are given a sequence ([Tj., k £ N) of random times on (U, A), and 
a sequence of random times (R^, k E N) on (S, C), which are such that for every k, 
Rk = Uk°Y holds true. Set Uoo = Roo = +oo, and consider the families of random 
variables 

y = {Y tARk , R u t e R+, k,le Noo}, X = {X tAUk , U h teR+,k,le N^}, 

where Noo = N U {+oo}. Then we obtain the following trivial statement, which for 
ease of later reference we record here as a lemma. 

Lemma 5.1. For every x E [0, 1], the family y has the same finite dimensional 
distributions under Q x , as X has under P x , that is, for all m, n £ N, t\, 
t m £ R+, h, i m , h, j n £ B h B m € S([0, 1] U {A}), d, 
...,C n £ S(M+U{+oo}), 

(5 5) ^ V** l/VRi l G Bl ' ' ' ' ' Yt ™ AR im G Bm > R h G Cl, . . . , Rj n £ C n ) 

= Px{X tl Mj n G Bi, . . . , X tmAUim £ B m , XJj x £ Ci, . . . , U jn £ C n ) 
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holds true. 

We shall denote by T x = (J 7 *, t £ R + ) the natural filtration of X, and set 

J* = a(X t , t G 

Below it will be convenient to work with both versions, X and Y, of the stochastic 
process simultaneously. The reason is that Y, having been constructed on product 
spaces, inherits a strong independence structure. This will turn out to be very con- 
venient to employ. On the other hand, on the path space SI = Ca(K+) we have the 
canonical families 6 = {Q t , t G 7 = (74, t G R+) of shift and stop operators, 
respectively: 

0t(w)(s) = w(s + t), s, t & M + , w S 
7t(u))(s) = w(s A i), s, t G M+, weQ, 

and we get 

X s o 9 t = X s+t , 

X s °lt = x sM . 

This implies that on Ca(K + ) we can use Galmarino's theorem (cf., e.g., flj p. 458], 
fl9] p. 86], [21] p. 43 ff], p4] p. 45]), which in particular gives a very useful charac- 
terization of the cr-algebra T x in the case when T is a stopping time with respect to 

T x : 

Lemma 5.2 (Galmarino's Theorem). 

(a) Let T be an M. + -valued random variable on (fi, ^4). Then T is an J rX - 
stopping time if and only if the following statement holds true: for all UJ\, 
W2 G O, t G R+, 7t(wi) = 74(^2)1 andT(uj\) < t imply T{oj\) = T(u2). 

(b) Let T be an T -stopping time, and let T x be the o-algebra of the T-past. 
Then T% = a(X tAT , t G M+). 

In order not to burden our notation too much, for both stochastic processes, Y and 
X, we shall denote the lifetimes by £ and the hitting times of sets C C M. by i^cs 
and if C = {x}, x G M, by for simplicity. Also, in both cases we shall write 
for the expectation E x {-), x G [0,1]. Whenever necessary, we add a corresponding 
superscript to avoid confusion. 

We bring in the notation 

c * fl, ifc = 0, 
° [0, ifc = l, 

and set T = 0, T x = H A For k G N, we define on the event {T fe < +00} and 
with c = Xx k 

{H c , if fc is odd, 
H c * o 9 Tk ; if A; is even. 

On {r fc = +00} put T fe+ i = +00. Obviously, S*. = T k o Y, k G N , and it 
follows that under Q x , x G [0, 1], Sk has the same law as under P x . Hence a.s. 
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the sequence (Tf., k G No) is strictly increasing to +00. Moreover, it is easy to 



check — for example with the help of lemma 5.2 — that for every k G No, Tfc is an 
J^-stopping time. 

Consider the process X defined as X being stopped at the stopping time T\, i.e., 



X with absorption in the endpoints 0, 1 of the interval [0, 1]. By lemma 5.1 we find 
that X is equivalent to Y with absorption in {0, 1}, and hence by the construction of 
Y, X is equivalent to a standard Brownian motion with absorption in {0, 1}. Then it 
is easy to check that X is a Markov process relative to T x . Futhermore, if instead 
of stopping we consider a process X which is defined as X with killing at the time 
Ti, we analogously get that X is equivalent to a standard Brownian motion which 
is killed at the boundary of the interval [0, 1], and X, too, has the Markov property 
relative to T x . 

5.2. Simple Markov Property. In this subsection we prove the following 
Theorem 5.3. X is a normal, homogeneous Markov process. 



The remainder of this subsection is devoted to the proof of theorem 5.3 We first 
prepare some lemmas, whose statements heuristically are quite obvious. Unfortu- 
nately, their formal proofs are rather technical, and therefore they are deferred to 
appendix [E] 

Throughout this subsection T = {Ft, t > 0) denotes the filtration F x of X, f is 
a bounded, measurable function on [0, 1] (extended to [0, 1] U {A} in the usual way 



via /(A) = 0), s, t belong to M + , and x to [0, 1]. The first two lemmas, lemma 5.4 



and lemma [53] below, form the key to the proof of theorem 5.3 



Lemma 5.4. X is a strong Markov process with respect to the stopping times T^, 
k G N. That is, for every bounded F^-measurable random variable V, all k G N, 
x G [0, 1], 

(5.7) E x {V o 9 Tk I F Tk ) = E XTk (V) , P x -a.s. on {X Tk + A} . 

We write X c for the process X with killing at c = or c = 1. Thus 

E x (f{X c t )) = E x (f(X t ); t<H c ), c = 0, 1. 

Similarly, A and B denote the processes A, B, with killing at 1, respectively. 

Lemma 5.5. Under P x , x G [0, 1], X 1 is equivalent to A under \i x , and X° is 
equivalent to B under fj, x . In particular, X 1 and X° are Markov processes. 

In appendix |E] we derive from these two lemmas the following three results. 

Lemma 5.6. On the event {0 < t < Tf, X t / A}, 

(5.8) E x (f(X s+t ) I Ft) = E Xt (f(X s )), P x -a.s., x G [0, 1]. 
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Lemma 5.7. On the event {X Tl = c, T k < t < T k+1 , X t ^ A}, k G N, the 
following holds true P x -a.s., x G [0, 1]: 

E x (f(X s+t ); T k <s + t<T k+1 \ F t ) 

= E Xt (f(X s );s<T 1 ) 

°' 9 ' 'E Xt (f(X s ); T 1 <s< T 2 , X Tl = c) , k odd, 



+ 



E Xt (f(X s ); T x < s < T 2 , X Tl =c*), k 



even. 



Lemma 5.8. On the event {X Tl = c, T k < t < T k+ i, X t / A}, k G N, the 
following holds true P x -a.s., x G [0, 1]: 

E x (f{X s+t ); T k+1 < s + t\T t ) 

(5. 10) _ J E Xt {f(X s ); T 2 < s, X Tl = c, or Ty < s, X Tl =c*), k odd, 

[E Xt (f{X s ); T 2 < s, X Tl = c*, orT x < s, X Tl = c) , k even. 



Now we are ready to prove theorem 5.3 Write 

oo 

E x (f(X s+t ) | F t ) =Y,Mf( x s+t); T k <t< T k+1 | Ft) 

k=0 

oo 

= ^E x {f{X s+t ) | F t ) l{T k <t<T k+1 }, 

k=0 

and therefore it is sufficient to prove the formula 

(5.11) E x (f(X s+t ) | Ft) = E Xt (f(X s )), P x - a .s. 

on each of the events {T k < t < T k+ \, X t / A}, k G No- For k = this is being 
taken care of by lemma [576] Hence we consider k G N from now on. On each of the 
events A fc = {T k < t < T k+ {\ we have T\ < t. This implies F Tl nA k cFtf] A k , 



and that Xt x is J^-measurable. This is so, because it is Ft x -measurable (e.g., 1 34 



Proposition 1.4.9], which can be applied because X being ^-adapted and having right 
continuous paths is progressively measurable relative to F). Therefore we get on 

{T k < t < Tfc+i, X t ^ A} 

E x (f(X s+t )\F t ) 

= E x (f(X s+t ); X Tl =0\ Ft) + E x (f(X s+t ); X Tl = l\ F t ) 



= E x (f(X s+t ) | Ft) l{x Tl =o} + Ex(f(X s +t) | F t ) l{x Tl =i}- 

Hence it is sufficient to prove equality ( |5.11| ) on each of the events {T k < t < 
T k+1 , X t + A} n {X Tl = 0}, and {T k < t < T k+1 , X t + A} n {X Tl = 1}, 
k G N. We shall only consider the first of these under the additional assumption that 
k is odd. The case when k is even, and the second term above can be treated in the 
same manner. 
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Fix k G N, k odd. On {T k < t < T k+1 , X t ^ A} n {X Tl = 0} we write 
E x (f(X s+t ) | T t ) = E x (f(X s+t ); T k <s + t<T k+l \ T t ) 

+ E x (f(X s+t ); T k+1 < s + t\ T t ) 



To the first of the terms on the right hand side we apply lemma 5.7 to the second 
lemma |5T8] and we obtain (P x -a.s.) 

E x (f(X s+t )\T t ) 

= E Xt (f(X s ); {T 1 <s< T 2 , X Tl = 0} U {s < T x }) 

+ E Xt (f(X s ); {T 2 < s, X Tl = 0} W {T x < s, X Tl = 1}) 

= E Xt (f(X s ); {s < Tt} W {Tx < s, X Tl = 0} « {T t < s, X Tl = 1}) 

= E Xt {f(X s ); {a<Ti}a{r 1 <fl}) 



which concludes the proof of theorem 5.3 



5.3. Strong Markov Property and Generator. In this section we prove the strong 
Markov property of X, and calculate its generator. 

Lemma 5.9. X is a Feller process. 

Proof. It is well-known, that it is sufficient to prove (i) that the resolvent of X pre- 
serves C([0, 1]), and (ii) that for all / G C([0, 1]), x G [0, 1], E x (f(X t )) converges 
to f(x) as t decreases to 0. (For example, a complete proof can be found in (24).) 
Statement (ii) immediately follows by an application of the dominated convergence 
theorem and the fact that X is a normal process with right continuous paths. 

Suppose that / G C([0, 1]), and consider the resolvent R = (R\, A > 0) of X. 
Because X is a strong Markov process with respect to the stopping times T k , k G N, 



(cf. lemma [5T4| > we get the first passage time formula (see, e.g., appendix |D|) for T\ in 
the form 

R x f(x) = E x (p e- xt f(X t ) dt) + E x (e- XT i R x f(X Tl )) 

,7i 



E X [J e- xt f(X t )dr 
+ E x (e- XTl ; X Tl = 0) R x f(0) + E x (e~ XT ^; X Tl = l) R x f(l). 



Lemma 5.5 states that X stopped at T\ is equivalent to a standard Brownian motion 
B which is stopped when reaching any of the boundary points of the interval [0,1]. 
Therefore we get for x G (0, 1) 

Rxf(x) = R^f(x) + E x {e~ XH o; ^ < Hf) R x f(0) 

+ E x {e- XH ?;H*<H*)R x f{l), 

where Hf, c = 0, 1, denotes the hitting time of c by B, and R D denotes the resolvent 
of B with Dirichlet boundary conditions at 0, 1 (which is equal to the resolvent of 
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expressions above are known explicitly (see, e.g., ||7||l9l): RP has the integral kernel 



X with Dirichlet boundary conditions at 0, 1, so that the notation is consistent). All 

2? 

■P( x ,y) = 1 ^2^ c -V2X\x-y+2k\ _ e -yM\x+v+-2l. ' 



2X kez 



while 



sinh(\/2A(l - x)) 



sinh(v / 2A) 



E x (e-™°;H? < H?) 



p ( P -\H?. h b , ttB\ _ sinh(y / 2Ax) 
E x (e i, Hl <H )- gmh(V _ ) 



It is now obvious that R\ maps C([0, 1]) into itself, and the proof is finished. □ 

Now we can apply standard results (cf., e.g, (34} Theorem III.3.1], or [39, Theo- 
rem III. 15.3]) to obtain 

Corollary 5.10. X is a strong Markov process relative to its natural filtration. In 



particular, X is a Brownian motion on [0, 1] in the sense of definition 1.1 



Remark 5.11. As discussed in, e.g., p0)[34][39) , on the basis of corollary 5.10 (and 
the path properties of X) one can prove that X is strongly Markovian relative to the 
canonical right continuous augmentation of its natural filtration. 

It remains to calculate the domain of the generator of X. We only discuss the 



boundary point 0, the arguments for 1 are similar. From lemma 5.5 we know that X 
with killing at 1 is equivalent to A with killing at 1. Therefore, starting in 0, if A has 
as a trap, then so has X, and if A jumps from to A after an exponential holding 
time, then the same is true for X. The remaining case is the one, where A and hence 
X leave immediately (into (0, 1)), and in this case we can use Dynkin's formula to 



compute the generator of X, because it is strongly Markovian (corollary 5.10). But 
this only involves the behavior of X before leaving an arbitrarily small neighborhood 
of 0, and hence at the generator of X is computed in the same way as for X with 
killing at 1, that is for the process A. Thus we have proved 



Theorem 5.12. X is a Brownian motion on [0, 1] in the sense of definition 1.1 and 
its generator is the second derivative on C 2 ([0, 1]) with boundary conditions at and 
1 given by ( |1.1[ ) and ( |1.2| ) respectively. 

Appendix A. Killing With A PCHAF 

This appendix provides a short account on the killing of a Markov process X = 
(Xt, t > 0) via a perfect continuous homogeneous additive functional (PCHAF) 
A = (A t , t > 0) of the process. Killing via the Feynman-Kac functional associated 
with A has been treated in detail for example in the books by Blumenthal-GetoorpJ 
and Dynkin (6j. Here we use a mechanism which has also been employed in |20 21| , 



and which is slightly different — in a sense somewhat more explicit than the method 
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of the Feynman-Kac functional just mentioned. In particular, we shall show how 
the simple and the strong Markov properties are preserved under this way of killing. 
This has also been argued in (2TJ for the case of Brownian motions, but we find 
the arguments given there quite difficult to follow. This appendix is based on the 
treatment in (3j of killing with the Feynman-Kac weight which, however, needs a 
number of non-trivial modifications. 

Throughout this appendix we shall use the terminology and notations of [3], except 
where otherwise indicated. (However, in order to be consistent with the other parts of 
this paper, there will be a few trivial variations in the notation, such as changing some 
superscripts into subscripts etc.) 

Assume that X = (Q,M,Mt, X t ,9t,P x ) is a temporally homogeneous Markov 
process with state space (E, £) as in section 1.3 of |3J. It will be convenient — and 
without loss of generality for the purposes of this paper — to assume that for all 
x 6 E, X has i-^-a.s. infinite lifetime. The natural filtration of X is denoted by 
F° = {F® , t > 0), its universal augmentation is denoted by F = (Ft, t > 0). In 
particular we have for all t > 0, F? C M t . We set F^ = cr(F^, t £ R + ), and 
similarly for Foo and .Moo- 



For the definition of a PCHAF we take the point of view of Williams [39 1, which 
has the advantage that the meaning of the word "perfect" is valid in the sense of 
both, Dynkin [6] (i.e., the additive functional is adapted to X), and of Blumenthal- 
Getoor (3j (i.e., the exceptional set where the additivity relation does not hold does 
not depend on the time variables): 

Definition A.l. A perfect continuous homogeneous additive functional (PCHAF) of 
X is an F- adapted process A with values in M + such that for some set Qq G Ai with 
Px(&o) = 1 far all x £ E, the following two properties hold for all u 6 Qq- 

(i) 1 1— > At(ui) is continuous from M + into itself non-decreasing, and Aq(u) = 0; 

(ii) for all s, t > 0, A s+t (u) = A s (uj) + A t o s (w). 

There will be no loss of generality for the sequel to assume that 1\) = Q: for 
example, for the purposes of this note one could put A t = 0, for all t > 0, on C^o- 
Moreover, we make the convention that for all A +00 (lu) = +oo. 

Remark A.2. Since t \— > At(uj) is continuous on R + for all uj € Q, we find that 
A is (8) v4/^(IR+)-measurable. Therefore, if T is any random time, i.e., a 

random variable with values in M + , then At{uj) = A T r w \ (uj) is a well defined random 
variable. Moreover, we also get the additivity relation in the form 

(A.l) A t+T ( u) (u)=A T fr ) (u) + A t o6 T ( u) (u), to en, ieIR + . 

Examples A3. Consider a standard Brownian family B = (B t , t > 0) on M. d , d £ N, 
and suppose that B has exclusively continuous paths. For example, one could realize 
B as the canonical coordinate mapping on Wiener space. Assume that V is any non- 
negative, measurable function on M. d . Consider 

(A.2) A t = [ V{B S ) ds, t > 0. 
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Then A is obviously a PCHAF. Moreover, in the literature (e.g., 1 19-21 34|) it is 
proved that for d = 1 and every x G M, the local time L t (x), t > 0, of B in x is a 
PCHAF. 

A.l. Killing. Define a new family of filtered probability spaces as follows. 

(A.3) n = n x !+, 

and denote a typical element in Cl by uj = (uj,s), w G (1, s £ M + . Define the 
following ex-algebras on Cl: 

(A.4) M=M®B{R+), 
and 

(A.5) A4° = .MxR + . 

In order to avoid that our notation becomes too unwieldy, we continue to denote the 
canonical extension of any random variable Z defined on (SI, M) to (SI, M) by the 
same symbol, i.e., we simply set Z(CS) = Z(uj), for uj = (oj,s) G 17. Clearly, for 
every t > the extensions of X t and A t are in A^°/£: and M°/B(R + ), respectively. 
For x 6 E set 

(A.6) P x = P x ®P e , 

where P e is the exponential law on i.e., it has density exp(— s) (s), 

s € E+. Furthermore for t > 0, cj € f2, a) = (w, s), define 

(A.7) 4 (w) = (0 t (uj) , s - A t (w) A s) . 

Then we have for all s, t > 0, X s o 6 t = X s+ t, and A s+ t = At + A s o ^. Moreover, 
if we set 

(A.8) S(u) = s, uj = (uj,s)en 
then for all i > 

(A.9) S o $ t = S - A t A 5. 
Define 

(A.10) C = inf{t > 0, At > 5}. 
Set 



(A. 11) X t 



\X t , if* < C. 
[A, ift>C, 

where again A is a cemetery point. 

Finally, for every t > we define a family 7W( C M. of subsets of O as follows. 
A set A G belongs to A^t if and only if there exists a set A t G A4 t so that 

(A.12) An{C>t} = (A ( xl + )n{C>i}. 

If A = {£ > t], t > 0, we only need to choose A t = f2 G A^ t to see that {( > i} G 
A4t, i.e., by construction of (Ait, t > 0), ( is a stopping time relative to (A4j , i > 0). 
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It is straightforward to show the following. 



Lemma A.4. (Ait, t > 0) is a filtration of sub— a— algebras ofAi, and X is adapted 
to this filtration. Furthermore, if(Mt, t > 0) is right continuous then so is (Mt, t > 
0). 

In the sequel we suppose that T is a random time defined on (Q,M), which as 
above is extended to (Cl, M.) in the canonical way. Let us define 

(A.13) {5> A T+ ) = \J{S> A T+e }. 

The proof of the following lemma is routine and therefore omitted. 



Lemma A.5. = S on the event {( < +oo}. Moreover, the following relations 
hold true: 



(A. 14) {C > T} = {S > A T }, 

(A.15) {(>T} = {S>A T+ }n{T<+oo}. 

Also the proof of the next lemma is left to the reader. 

Lemma A.6. Let T be a random time as above. Then on {Q> T} 

(A.16) ( = T + (o§ T . 

Lemma A.7. Let T be a random time as above, and let x £ E. Then 

(A. 17) P X (C > T\M°) = e~ Ar 

holds true P r -a.s. 



Proof. Observe that the right hand side of equation ( A.17| ) is A4° measurable. Let 
A G .M , i.e., A = A x M + with A G M. Without loss of generality we may assume 



that A C {T < +oo}. Using Fubini's theorem, lemma A.5 and the continuity of 
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1 1 — y At{uj), to £ Q, on M + , we can calculate in the following way: 

P x ({( > T} n A) = P X ({A T+ <S}nk) 

= HmP B ({4 r +e<5}nA) 



- 1 [i4T+.,+oo]°<Sr4 p *® i «) 

40 JAxK + 



CO 



= lim / / e _fl ds dP x (w) 

= lim£; a: (e" j4T + e ;A) 
= ^(e- AT ;A) 

= E x (e- AT ;A). □ 
Recall our convention that every function / on E is extended to with /(A) = 

0. 

Corollary A.8 (Feynman-Kac-formula). For every bounded measurable function f 
on E, and all t > 0, x 6 i£, the following formula holds true 

(A.18) E x (f(X t )) =E x (f(X t )e~ A *). 



Proof Since /(A t ) G A4°/B(R), we get from lemma \AJ\ (with the choice T = t) 

E x (f(X t ))=E x (f(X t yX>t) 

= E x (f(X t )P x ((>t\M )) 
= E x (f(X t )e- At ) 

= E x (f(X t )e- A <). □ 

A. 2. Preservation of the Markov Property. Now we prove the 
Theorem A.9. X is a Markov process with respect to (Ait, t > 0). 



Proof. Let / be a bounded, measurable function on E. We remark that corollary A.8 
shows that x h-> E x (f(X t )) is £*/£>(M)-measurable, where £* denotes the universal 
augmentation of £. 

Let x £ E, s, t > 0. We have to prove that 

(A. 19) E x (f(X s+t ) | Mt) = Ex t (f(X a )), A-a.s. 
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To this end, let A e Mt, and choose A t G M t to be such that A fl {( > t} 
(At x M + ) n {( > t}. Then we can compute in the following way: 

E x (f(X s+t );A) = E x (f(X s+t ); An{(>s + f}) 

= E x (f(X s+t ); (A n {C > t}) n {C > s + t}) 

= E x (f(X s+t ); (A t xl + )n{(>s + t}) 



E x (f(X s+t ) l AtxS+ 
E x (f(X s+t ) l AtxS+ 
E x (f(X s+t )e- A °+t 



P x ((> s + t\M )) 

;A t ), 



where we made use of lemma A.7| By the additivity of A the last expression is equal 
to 

E x ((f(X s ) er A *) o t e~ Ai ;A t ) = E x [E x ({f{X s ) e~ A °) o 6 t \ M t ) e~ At -A t ) 

= E x (E Xt {f(X s )e- A °)e- A *;A t ) 
= E x (F(X t )e- At ;A t ), 

where we have used the Markov property of X relative to (Alt, t > 0), and where 
we have set 



(A.20) 



F(x)=E t (f(X,)e- A '). 
F(x) = E x {f(X s ). 



Observe that by corollary |A.8| 
(A.21) 
Thus we have 

E x (f(X s+t );A) =E x (F(X t )e- At ;A t ) 



E x {F(X t )l 



-At 



A t xl 



) 



EjF(X t )l AtxW+ P x {(>t\M 



E x [F(X t ] 
E x {F(X t ] 
E x (F(X t ] 



(A t x M+) n {( > t}) 

An{C>t}) 
A) 



4 (%,(/(£)); A), 



and the proof of relation (A.19> is complete. 



□ 



A.3. Preservation of the Strong Markov Property. In this subsection we consider 
the strong Markov property. So we assume in addition that X is a strong Markov 
process as in section 1.8 of j3j. 

The crucial point is the following analogue of Lemma III.3.9, p. 108, in (3j. 
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Lemma A. 10. Let T be a stopping time relative to (Ait, t > 0). 

(a) There is a unique stopping time T relative to (Ait, t > 0) so tha^ 
(A.22) f A C = T A C 

holds. 

(b) Let T and T be as in (a), and let A G Aif. Then there is a set A G Aix so 
that 

(A.23) An{C>T} = (AxI + )n{C>T}. 

Proof. The proof of statement (a) is very similar to the corresponding one in Lem- 
ma III.3.9 in j3j, and is therefore omitted. 

We begin the proof of statement (b) with the following remark which will prove to 
be useful. Let r > and w £ fibe given. We claim there exists ui = (uj, s) £ Cl, 
s G IR + , so that ((uj) > r. In fact, we note that by equation ( |A. 15| ) this claim 
is equivalent to saying that there is an e > so that A T+e (co) < S(uj). Choose 
any e > 0, and observe that the continuity of the increasing mapping u \-t A u (uj) 
guarantees the existence of s > so that s > A T+e (u). Set to = (u, s). Then the 
claim follows because S(u)) = s. 

Now let T be a stopping time with respect to (Aif, t > 0), and let T be the 
associated (Ait, t > 0)-stopping time as in part (a). Choose A G Aif, and for t > 
define the set 

A t = An{f < t, C > *}. 

Since A n {f < t} G M t , we have that 

(A.24) A ( = (A i xl + )n{(>t}, 

for some set At G Ai t . We show that A f C {T < t}: If uj G A t pick s G M+ so that 
((u, s) > t (which is always possible due the remark above). Therefore 

(u,s) G A t C {f < t, ( > t}. 



Equation ( |A.22| ) entails that 

(A.25) {f < t, C > t} = ({T <t}x 1+) n {C > t}, 

and therefore we get (uj, s) G {T < t}. 

Similarly one shows that for t', t > with t' < t, one has Af C At, where these 
two sets are associated with At and A t / as in ( A.24 1. 

Now we claim that for t', t > with t' < t 

Af n {T < t'} =A t n{T< t'} 
holds true. Since we have that Af C At, we only need to show the inclusion 

Af n {T < t'} d A t n {T < t'}. 



7 

Recall that we denote the canonical extension of T from Q to £1 again by T. 
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To this end, assume that oj G Af n {T < if}. Choose s > so that (u, s) G {C > 
t} C {C > f }. Then 

(w, s) G (A t x 1+) n {T < t} n {C > t) = A t . 



Therefore (u, s) G A n {T < t' , ( > t'}, i.e., (w, s) G A t /. Using formula ( [Agg ) for 
A^, we get that ui G Af Pi {T < t'}, and our claim is proved. 

We conclude that for every t > 0, At G .Mt- Indeed, let < t, then 

A* n {T < t'} = At/ n {T < t'} = Af G A^f . 

If t' > t, then A t C {T < t} C {T < t'} and therefore 

A t n{T <t'} = A t £ M t C M t >. 

To finish the proof of the lemma, set 

A = |J A u 

where Q + = Qfl M+. Then we find that A G M-t, and moreover by what has just 
been proved 

An{T<t} = A t , t G Q+. 
Now we can calculate as follows 

A n {C > T} = A n ( (J {f < t, C > t}' 
= |J An{f <i, c>t} 

t6Q+ 

= |J (A 4 x R+) n{(>t} 

= \J (A x M + ) n ({t < t} x 1+) n {C > t} 

= (A x R + ) n ( |J {f < t, c > *}) 

= (AxR + )n{C>t}, 
and the proof of the lemma is complete. □ 

For the statement of the strong Markov property of X, we need the following 
lemma whose proof is straightforward in view of lemma A. 10 and therefore left to the 
reader: 

Lemma A.ll. Let T be stopping time relative to (A4t, t > 0). Then Xf G M.f/£. 
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Theorem A.12. Suppose that X is a strong Markov process relative to (Ait, t > 0). 
Then X is a strong Markov process relative to (Ait, t > 0). 

Proof. Let T be an (Ait, t > 0)-stoppin g time , and let T be the associated stopping 
with respect to (Ait, t > 0) as in lemma A. 10 Let x G E, A G Aif and consider a 



bounded measurable function / on E. According to lemma A. 10 b there exists a set 
A G Air so that relation ( |A.23[ ) holds true. With a similar calculation as in the proof 
of theorem A.9| and lemma [A77] we get 

E x (f(X s+f );A) = E x (f(X s+T )e- A °+ T ;A), s > 0. 
With equation ( |A.1[ ) and the assumption that AT is a strong Markov process we obtain 

4 {f(X s+ f); A) = E x ((f(X.) e~ A *) o T e -^ ; a) 

= ^(i? XT (/(X s )e-^)e-^;A) 

= i(%(/(I s )e^)e^;Ax 
With equation ( |A.17[ ) and the notation ( |A.20] > we find 

E x {f(X s+f ); A) = 4 (f(X t ) P x (C>T\M°);Ax 
= E x (f(X t );(A xR + ) n {(> T} 



where we used the fact that the indicator of A x M + and F(Xt) are both Ai°- 
measurable. Because of f = T, Xf = X T on {( > T} = {( > f}, and due 
to relation ( |A.23| > we get 

E x (f(X s+f );A) =E x (F(X f )-A) 

= E x (E jtf (f(X s ));A). 

As a consequence, P^-a.s. 

(A.26) E x (f(X g+f ) | M f ) = E^_ (f(X s )) , 

which concludes the proof. □ 

Appendix B. Some Results related to Brownian Local Time 

We consider a standard one-dimensional Brownian motion B = (B t , t > 0) 
on a family (Q,A, (P x , x G M)) of probability spaces, endowed with a right con- 
tinuous, complete filtration F = (Ft, t > 0). For given x G M we denote by 
L(x) = (L t (x), t > 0) its local time at x. We assume, as we may, that both have 
exclusively continuous paths. (For example, for the local time L(x) this follows from 



the Tanaka formula (201 Proposition 3.6.8] for L(x), and Skorokhod's Lemma |20 



Lemma 3.6.14].) In particular, for given x G M, L(x) is a progressively measurable 
stochastic process. The right continuous pseudo-inverse of L will be denoted by K: 

(B.l) K r = inf {t >0,L t > r}, r G R+, 
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with the convention inf = +00. Because L is ^-adapted and F is right continuous, 
for every r > 0, K r is an ^-stopping time. 

Lemma B.l. For r > 0, 

(B.2) P (K r G dl) = e- r ' 2/21 dl, I > 0, 

Y^27r/ 3 

and 

(B.3) ^ (e" AKr ) = e"^ r , A > 0. 



Proof. By Tanaka's formula (e.g., | |20] Chapter 3.6]), for every q > the stochastic 
process 

t ^ e -a{\Bt\-L t )-a 2 t/2 

is a martingale under Pq. Since L t f +00, P - a - s -> we g et -Fb(-Kr < +00) = 1. The 
usual argument of truncation of K r , application of Doob's optional sampling/stopping 



theorem (e.g., |20, p. 19], or |34J p. 65, ff.]), followed by the removal of the truncation 
yields with a = y/2X, 

E (e- XKr )=e W ** r , r>0, A > 0. 
Formula ( |B.2| ) follows from equation ( |B.3| > by inversion of the Laplace transform with 



the help of the formula ( |3.15| >, and the uniqueness theorem for Laplace transforms. 

□ 

Let a > 0, and consider the stopping time H a , defined as the hitting time of a by 
the reflecting Brownian motion \B\, i.e., H a is the hitting time of the set {—a, a} by 
B. We want to calculate the law of Ln a (0) under Pq with a method similar to the 
proof of Theorem 6.4.3 in [20] — actually we shall solve a variant of Problem 6.4.4 
in (20|. 

Denote by C 2 (IR; 0, ±a) the space of bounded, continuous functions / on K which 
belong to C 2 (M \ {—a, 0, a}), and which are such that at —a, 0, a, f and /" have 
(finite) left and right limits. For / G C 2 (M; 0, ±a), t > 0, the generalized Ito formula 
reads (cf. [20, Equation 3.6.53]) 

df(B t ) = f(B t )dB t + \f"{B t )dt 

(B.4) 

+ - (A/'(0) dLt(0) + Af'(-a) dL t (-a) + Af'(a) dL t (a] 

where we denoted Af'(b) = f'(b+) - f(b-), b G R. For a, f3 > 0, p G E set 

(B.5) Z t = exp(-at-^ t (0) + p(Lt(-a) + Lt(a))), 

and note that for all t > 0, Z t belongs to all C P {P X ), p > 1, a; G K. (For example, 
this follows easily from an application of Holder's inequality and of the strong Markov 
property of B, the translation invariance of the transition density of B, and the explicit 



law of L t (0) under Pq, e.g., Theorem 3.6.17 in [20 1.) 
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Now let v G C 2 (M; 0, ±a) be an even function, and define the stochastic process 



(B.6) 



X t = v{B t )Z u t>0. 



We choose v in such a way that X becomes for P x , x G K, a martingale, namely, we 
require that u solves the following boundary value problem: 

1 „ 

-v = av, 

v'(0+) = (3v(0), 
v(a) = 1. 



(B.7a) 

(B.7b) 
(B.7c) 



on 



\ {-a,0,a}, 



Furthermore, in the definition ( |B.5 1 of Z we choose 

1 



(B.8) 



Av'(a) 



Then it is straightforward to check with formula ( |B.4 1 that Xt is indeed a martingale 
under P x , x G R. 

We remark that L t (±a) = for all £ < i7 a . Now we proceed as in the proof of 
lemma B.l We truncate H a so that it becomes a finite stopping time, apply Doob's 
optional sampling theorem, and remove the truncation by an application of the domi- 
nated convergence theorem. The result is the formula 



(B.9) 

and in particular 
(B.10) 



i{x) = E x (e 



- a H a -(3L Ha (0) 



x G 



E {e 



limf (0). 

c40 



(B.ll) 



v(x) 



The unique bounded even solution of the boundary problem ( |B.7[ ) is given by 

2a cosh(v / 2ax) + (3 sinh.(^/2ax) 
2a cosh(\/2cra) + P sinh(\/2cra) 

2a(x— a) 



< X < a, 
x > a, 

1 



and v(x) = v(—x), x < 0. Thus we get 
(B.12) E (e-P L »*W) 
and we have proved the 



Lemma B.2. Under Pq, Ljj a (0) is exponentially distributed with mean a. 

Remark B.3. For the reader who compares this result with equation (6.4.16) in [20], 
we want point out that here H a denotes the hitting time of {—a, a} while there T& is 
the hitting time of b, and furthermore our convention for the local time differs by the 
one used in 1 20 1 by a factor 2. 
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For the next result we continue to denote by H a the hitting time of the set {—a, a} 
by the Brownian motion B, and use the set-up of appendix [A] Fix f3 > 0, and denote 
by P/3 the exponential law on (R, 0(R)) of rate /3. Construct the family of product 
spaces (Cl, A, (P x , x £ M)) analogously as in appendix |a| with P e replaced by Pp 
in equation ( |A.6| ). For a PCHAF we choose here the local time L(0) at zero, and 
consider the stopping time Qp defined as in ( |A.10[ ) with A t = L t (0), t > 0. 

Lemma B.4. For all a, (3 > 0, 

(B-13) P Q {H a < Qp) = —^-j- 

' 1 + pa 

Proof. We use equation ( A.14 1, and lemma B.2| to compute 

Po[H a < (p) = Po[S > L Ha (0)) 

= (3 e-? s P {L Ha (0) < s) ds 
Jo 



ds 



1 



1 + 13a 

The proof is concluded by the remark that by definition of H a and C we have that 
H a ^ £p unless H a = C/3 = +oo, but the latter event has Po _measure zero. □ 

Appendix C. An Inverse Laplace Transform and Heat Kernels 

The following lemma gives a generalization of a well-known formula for Laplace 
transforms (e.g., (51 eq. (5.1.2)]. 

Lemma C.l. Suppose that f is a bounded, measurable function on M+, and let a > 0, 
b > 0. Define 

(C.l) K a uf{t) = .f± I* * + 6 q/ , expf-a f( s ) ds. 



2tt Jo (t-sf/ 2 *A 2(t-s) 
Then 

(C.2) C(K a<b f)(X) = e- V2 ^ b Cf(X + V2a~X), A > 0, 

where C denotes the Laplace transform. 



Proof. Let A > 0. The theorem of Fubini-Tonelli yields 

Va (s + b) 

V27r(i - s) 3 ^"V ~2(i-s) 



C(K a , b f)(X)= I /(.) /°V* W + h \ exp(-afe^W 



The inner integral is equal to 

(C 3) n e~ X{t+s) V^(l±j) e -a(s+6) 2 /2* dt = e ~\s e -V2^\(s+b) 

Jo V2tt^ 
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where we used the Laplace transform (|3.15|>. Inserting the right hand side of equa 



tion ( |C.3| > above, we find formula ( |C.2[ ). □ 



Recall the definiti on (|4~6] > for p(X), A > 0, and ( |4.14| ) for gp n (t,x), /3, 7, t > 0, 



x > 0. Then lemma C.l immediately gives with the choice a = 7 , b = jx, and 



f(s) = exp(— ^s/7) the following result. 

Corollary C.2. For /3, 7, x > 0, f/ze inverse Laplace transform of 

(C.4) A i — ^ p(A) e - ^, A > 0, 

« g/ve« by t gp t j(t, x), t > 0. 

In order to relate 7 to other heat kernels in this article we consider the limits 
when the parameters (3, 7 tend to zero. First fix (3 > 0. Then a straightforward 
(though somewhat tedious) calculation shows that for all t > 0, x > 0, 

(C.5) lim ag 7 (i, x) = ag (t, x), 

74-0 



where gpfi(t, x) is given by ( |2.24 ). On the other hand, for fixed 7 > 0, one has 



(C.6) Km%(t,!i;) =50,7(^,2;), 

with go^(t, x) given by ( |3.23| ). Maybe the easiest way to see this is to compute the 



limit of the expression defining gp~, i.e., formula ( |4.14| ), and then to apply lemma C.l 
to show that the Laplace transform of the resulting integral is given by 

A^^ e"^, A>0. 

^/2A + 7 A 

Now one can invert this Laplace transform, e.g., with formula (5.6.16) in [8], and this 
gives go,7 as defined above. Finally, with formulae ( |2.24[ ), ( |3.23| ) it is easy to see that 

(C.7) )wigp Q {t,x) =YangQ tl (t,x) =g(t,x). 

P4.0 ■yiO 

Appendix D. First Passage Time Formula 

For the convenience of the reader we provide in this section a proof of a specific 
form of the well-known first passage time formula (e.g., [ [19| p. 94], (2T| p. 158]) 
which we use in this article. 

We assume that X is a normal, strong Markov process relative to a right continuous 
filtration T, with state space (E, £) being a locally compact, separable metric space, 
such that X has cadlag paths and F is complete with respect to the family P = 
(P x , x G E) of probability measures on the underlying measurable space (Q, A). 

Let xq be a given point in E, and denote by Ho its hitting time by X. We suppose 
that P-a.s. Hq is finite. In particular, this implies that X cannot be killed before 
reaching xq. Moreover, our assumptions above entail that Hq is an ^-stopping time 
(cf., e.g., Theorem III.2.17 in [34]). The resolvent of X will be denoted by R = 
(R\, A > 0), while R° = A > 0) stands for the resolvent of the process X° 
which is obtained from X by killing it when reaching the point xq. Recall that by 



46 V. KOSTRYKIN, J. POTTHOFF, AND R. SCHRADER 

convention every real valued function / on E is extended to E U {A} by setting 
/(A) = 0. Then we have the 

Lemma D.l. For every bounded measurable functions f on E, and all x G E, A > 0, 
(D.l) R x f(x) = R° x f(x) + E x (e- XHo ) R x f(x ). 

Proof. Let /, A and x be as in the hypothesis. Then 

R x f(x) = E X (J e- xt f(X t ) dt) 

Ho \ / f<x> 

— \tf/V\i,\ . t-t I / -—At. 



t+H 



\dt 



= E x [j o e-^fiX^dtj+E^j^ e- M f(X t )dt 
= E x (^°° e~ xt f{X Q t ) dt) + E x (e~ XH ° e- xt f(X, 
= R°J(x) + E x (e- XH ° e- xt E x {f(X t+Ho ) \ F Ho ) dt 

= R\f(x) + E x (e~ XH ° e~ Xt E X0 (f(X t )) dt) 

= R°J(x) + E x (e- XH °) Rxf(xo), 

where we used the strong Markov property of X, and its right continuity to conclude 
that Xjj = xq. □ 

Appendix E. Proofs of Lemmas [5T4T - I5 . 81 
In this appendix we prove lemmas |5.4f|5.8| 

Proof of lemma \5A\ Fix k G N, x G [0, 1]. By standard arguments it is enough to 
show equality ( |5.7 ) for V of the form 

n 

V = Ufi(X ti ), 
1=1 

with n G N, < t\ < ti < ■ ■ ■ < t n , and bounded measurable functions fx, /2,. . . , 
/„ on [0,1]. Thus we have to show that 

n n 

(E.1) E x (j{f t {X u+Tk )- A)=E x (E XTk (j[f i (X ti )); a) 

i=\ * i=l 

for every A G H {XT k 7^ A}. It is sufficient to show ( |E.1[ ) for every A in a in- 
stable generator of J 7 ^ n {Xx k 7^ A}, and since by Galmarino's theorem, lemma 5.2 
J 7 * = a(X , AT k ), we can choose A to be of the form 



A=f]{x S]ATk eC j }n{x Tk ^A}, 
3=1 
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with m G N, si, s 2 , s m G R+, Ci, C 2 , . . . , C m G S([0, 1] U {A}). Remark 
that by construction {Xx k ^ A} = G {0, 1}}. Hence we may decompose 

A = A l±l Ai in T$ n {A Tfc / A} with A c = A n {A Tfc = c}, c = 0, 1. Consider 
the right hand side of equation ( |E.1| ) with A replaced by Ao- By lemma 5.1 we find 

n n 

(E.2) E x (jlf t (X ti+Tk )-Ao) = E x (jlf i (Y ti+ s h y,A% 

i=l ' i=l 



with 



^ =C]{Ys J AS k ^C J }n{Y Sk = 0}. 

3=1 

In order to compute the expectation value on the right hand side of equation (E.2), we 
factorize the probability space (S, C, Q x ) as follows: 



i-.<n-l 



n— 1 oo 

_ ^/ ^>n _ ^ -/ 

fc=0 fc=n 
n— 1 oo 



^ = 0^, C^ n = (g)C fc , 

fe=0 fc=n 

71—1 OO 

Q|»-l = QO ® ((g) , Q>™ = (g) 



fc=l fc=7l 



We show that A^ belongs to the er-algebra C-™" 1 x 3- n , i.e., that the indicator of 
this set only is a function of the variables in H- n_1 . Firstly, from the construction of 
Y and the sequence [Sk, G N) it follows that 



{Y Sk = 0} 



{Pi < Po) n {^fc < if A; is even, 

{Po < P?} n {S^ < +oo}, if A; is odd. 



Since in both cases the sets on the right hand side belong to C- n 1 x H- n , we find 
that {Y Sk = 0} G C^ n-1 x E^ n . Next, for j G { 1 , . . . , m} we get 

{Y Sj AS k eC 3 }n{Y Sk =0} 

= {y Sj g c,-, Sj < s k , Y Sk = o} a {y Sfc g c,-, ^ > s fc , r 5fe = o}. 

The first term belongs to C- n_1 x H- n , and for the second we note that it is empty 
unless G Cj, and in this case it is equal to {sj > Sk} H {Ys h = 0}. Therefore this 
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event belongs to C- n ~ x x H- n , too. Hence we can write 

n 

E x (j[MY tt+Sk );Al) 

= Ln-i 1a " OL n WW*) dQ^ n ) dQ^- 1 

(E-3) = E / V^-^-' 



fe<Zi <•»<{„ 



1<n-l 



with 



(E.4) i2 Jl ,..., In = / Il/iC^+sJl^^^}^ 

and we made use of the fact that Q^-a.s. the sequence (S^, n S N) is strictly increas- 
ing to infinity. 

We claim that on {Ys k =0} 

n 

( R5 ) R k,-,ln = E o(Y[fi(Y ti ) l{^_ w . 1 <t i <S',._ Jk+a }) 

1=1 

holds true. Consider the integrand on the right hand side of equation ( |E.4[ ), choose 
i € {1, . . . , n}, and assume first that l{ = k. Then (on {Yg k = 0}) the corresponding 
term in the product is equal to 

fi(Xu+S k ) l{0<U<a k } = fi{ A u) !{0<U<a k }- 

Note that the term on the right hand side is constant on {Ys k = 0} C H- n_1 , and that 
under Q- n it has the same law as 

under Qq. Next consider the case ij > k + 1. By construction the relations 

S k -S k = a k + T k+1 + --- + V li ~\ 
(hX}) S h+1 -S k = a k + r k+1 + --- + V h , 

hold true on {Y$ k = 0}, where for r G N, r > k, 

J <7 r , if r — /c is even, 
1 r r , if r — k is odd. 

Moreover, on {Y Sk = 0} n < + S fc < S h+1 }, 



(E.7) 



(E-8) Y ti+Sk 



^_ {Sl -S k y if^-^iseven, 
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Remark that the right hand sides of the relations ( |E.6| >, ( |E.8| > are constant on H 
On the other hand, under Qq, i.e., with start of Y at 0, we get pfj = 0, and hence 

S h - k +i = a 1 + r 2 + --- + 

S h - k+ 2=<T 1 +T 2 + ---+^- k+1 , 
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(E.9) 

with (r e N, r > 1) 
(E.10) 



|cr r , if r is odd, 
lr r , if r is even. 



Furthermore, if U G [^-fc+i, then 

. ^Zi-fe+i if L - A; is even, 



(E.11) 



5f i_ £ +1 , if U - k is odd. 



A comparison of the right hand sides of equations (E^MEJS) with (E.9 H E.ll ) 
shows that on {Ys k = 0} formula ( E.5 1 holds true as claimed. 
We insert formula ( |E.5| ) into equation ( |E.3| ) and get 



7t 

^(n/i(^+s fe );A 

1=1 



X] ( 5^( A o') 1 {s, i _ fc+ i<ti<s li _ fc+2 } 

k<h<-<l„ 



i=l 



n 

= Q x (A u ')^o(n^(^) 1 {5 1 < ti ; 

i=l 
n 

= Ox(A2 r )^(n/i(y ti : 

t=l 

n 

= ^(^ (n/i(^)); A r 



j=i 



n 

= ^(^y Sfc (n^(^));A 



where for the third equality we used that under Qo, Si = 0. With lemma 5.1 we find 

n n 

E x (J] fi&U+n); Ao) = (#x Tfc ([J /i(X ti )) ; Ao) . 
i=i ' i=i 

We add this equation to the corresponding one where Ao is interchanged with Ai, and 
obtain 

n n 

E x ([J fi(X ti+Tk ); A) = E x (E XTk (j[ /,(X ti )) ; A 

i=l " i=l 
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Thus lemma 5.4 is proved. □ 

For later use we bring the strong Markov property of the last lemma into an efficient 
form. To this end, let us denote the family of cylinder functions of X by Q, i.e., G G Q 
is a random variable of the form 

(E.12) G = g(X n ,...,X Sm ), 

where m G N, g is a bounded, measurable function on [0, l] m (extended to ([0, 1] U 
{A}) m in the usual way), s\, s m G M+, with s\ < S2 < ■ • • < s m . For G G Q of 
this form set 

(E.13) G r = g(X Sl - r , . . . ,X Sm - r ), T < S X . 

Lemma E.l. Assume that G G Q is of the form ( |E. 12 1. Suppose furthermore that 
k G N, and that (p is a bounded, measurable function on R + . Then for all x G [0, 1], 

E x (cp(T k+1 )G; T k < 8l \M) 

(E.14) 

= E x Tk [<p(r + H c ) G r ) l [0 , ai ](r) | r=TfejC=(XTfc) » 
holds true P x —a.s. on {Xx k ^ A}. 

Proof. Let A G Pi {Xr fc / A}, and decompose it in J 7 ^ as A = A W Aiwith 
A c = An {X Tfc =c},c = 0,1. 

We choose first </?(x) = exp(zax), a G M, a; G M, y?(+oo) = 1. Set uj = 
Sj — si > 0, j = 1, . . . , m, and define for s > 

-^( s ) = g{x s ,x U2+s , . . .,x Um+s ), 

so that G = K(si). We compute the following Laplace transform at A > 0: 



poo 

/ e- Xs E x (ip(T k+1 )K(s); A D {T k < s}) ds 
Jo 

POO 

= E x U aTk + 1 e- Xs K(s)ds; A 

J Ti. 



= E x (e iaT ^- XT « (^°° e- Xs K(s) ds} o Tfc ; A ). 

By construction T k+ \ = T k + H\ o 9r k when Xx k = 0. Therefore the last expression 
equals 

"OO 



E x ^ ia - x ^E x (^e iaH ^ j\- Xs K{s)ds) o9 Tk \ F Tk ); A ] 
= E x (e {ia - x)Tk - A ) E (e iaHl e" As K(s) ds} 

f'OQ f'OO 

= / e- Xs e ias P x {{T k G ds} n Ao) / e- As £ (e iQ ^ ds, 
Jo Jo 
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where we used the strong Markov property of X relative to Tk, cf. lemma 5.4 Invert- 
ing the Laplace transform we get 

E x (e iaT ^ G; A n{T k < Sl }) 



f Sl P X (A n {T k e dr}) E {e ta( ~ r+H ^ G r ) 
Jo 

/•CO 

/ P x (A n {T k e dr, T k < si}) E ( e ia ( r+H ^ G r ) 
Jo 

E x (E {e^+ H ^ G r ) l [0)Sl] (r) | r=T j A ) 
E*(Ex Tk (e**^ G r ) l [0)Sl] (r) | r=Tfe , c=(XTfc) ,; A 



Repeating the same calculation with Ao replaced by Ai, gives an analogous expres- 
sion, and adding both we obtain the following formula 

E x (e iaT *+i G; An {T k < si}) 

= E x (E XTk (e-(^) G r ) l [0>si] (r) \ r=ThtC={XT y, a) 

for arbitrary A G F x k n {Xr k ^ A}. Finally, a routine approximation argument 
combined with Fourier analysis yields equation ( |E.14| ). □ 



Now we are ready for the 



Proof of lemma\5^6\ We work on the event {0 < t < Ti, X t A}, t G R+. For 
s > 0, write 

E x (f(X s+t )\P x )=E x (f(X s+t ); s + t<Tx\F x ) 

(E-15) / 

+ E x (f{X s+t );s + t>T 1 \F x ). 

The first term on the right hand side only involves the process X with killing at the 
boundary of the interval [0, 1], and this process is Markovian relative to F x (cf. the 
discussion at the end of subsection |5.1[ ). Therefore this term is equal to 

(E.16) E Xt (f(X s );s<T 1 ). 

For the second term on the right hand side of equation ( |E.15[ ) we remark that 

F? n < Ti} c f£ n {t < Ti}, 

and therefore it is equal to 

E x (E x {f(X s+t ) | F%) ; Ti < s + t, X Tl = | F t x ) 

+ E x [E x {f(X s+t ) \F? 1 );T 1 <s + t,X Tl = l\F x 

= E x (MTi); Tx < s + t, X Tl = | F x ) 

+ E x {<px(T x ); T x < s + t, X Tl = 1 1 jf ), 
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where 

ip c (r) = E c (f(X s+t -r)) , c = 0, 1, < r < a + 1, 

and we used the strong Markov property of X relative to T\, lemma [EJ] The last 
two expectation values only involve the process X, that is, X with absorption at the 
endpoints of the interval [0,1], and its absorption time T\. T his i s a Markov process 
with respect to T x (cf. our discussion at the end of subsection 5.1 1. Moreover, t <T\ 
entails the relations cp c (Tx) = ip c (t + Tx) o t , {Tx < s + t} = {Tx o 6 t < s}, and 
{Xt x = c} = {Xt x o 9t = c}. Thus we get 

E x [f{X s+t );T x < S + t\FX) 

= E Xt (ip (t + Ti); Ti < s, X Tl = 0) 

+ E Xt ( Vl (t + Ti); T x < s, X Tl = 1) 

= E Xt (E (f(X s - r )) \ r=n , T x < s, X Tl = o) 

+ Ex^E^fiX^)) | r=Ti ; Tx < s, X Tl = l] 
= E Xt (Ex Tl [f(X s „ r ))\ r=Ti ;T 1 <s^ 
We apply lemma [ET] to the last expression, and the result is 

E Xt (E Xt (f(X s ) | Jg) ; ^ < s) = E Xt (f(X s ); T x < s) . 



Combining this with ( |E.16[ ) for the first term on the right hand side of equation ( E.15 1, 

we get on {t < Tx, X t / A} 

E x (f(X s+t ) | Tf ) = E Xt {f(X s )), P x - a .s., 



and lemma 5.6 is proved. □ 

Recall that we write X c for the process X with killing at c, c £ {0, 1}, and X 
for the process X with killing at either of the points 0, 1. Furthermore, A and B 
denote the processes A, B, with killing at 1, respectively. We shall make use of the 
semigroups U 1 , U , generated by A, B respectively: 

U}f{x) = E x [f(A t )) = E x (f(A t ); t < H?), 
U?f(x) = E x [f(B t )) = E x (f(B t ); t < F B ). 

With the help of the strong Markov property of the processes A, B, it is not hard 
to prove the following result: 

Lemma E.2. If s < t, then a.s.for all y € [0, 1], 

E y (f(A t ); t<H?\ Tf) = Ulj{A s ) 1 {s<h a } , 

(E.17) D o n 

E y {f(B t ); t<H*\ Jf ) = Ulj(B s ) 1 {s<h b } . 
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Proof of lemma 5.5 We only show the equivalence of X 1 and A, the equivalence of 
X° and B follows from an analogous argument. Also in this proof, for the sake of 
simplicity of notation we denote by F the filtration generated by A, and its hitting 
times of c = 0, 1, by H c . Fix t > 0, let G be a bounded J^-measurable random 
variable, < s < t, and suppose that F is a bounded measurable function on [0, 1]. 
Since {Hq < H{\ € Th (see, e.g., p0| Lemma 1.2.16]), we get with the strong 
Markov property of A, 

E x {GF(A t y,8<H <t<H 1 ) 

= E X [G 1{ s< h <h 1 } F(A t ) l{t<Hi}) 

= E X (G l {s<Ho<Hl y E x (F(A t ) l{ t< Hx} \ Fh )) 

= E x (Gl {s<Ho<Hl} E (F(A t „ r ); t-r<H x ) l [0)t] (r) \ r=H J, 

where we can, for example, use the Laplace transform as in the proof of lemma [Err] 
together with fact that on {Hq < Hi}, we have H\ = Hq-\- H\o 6 Ho , to prove the 
last equation. Therefore we obtain the following formula 

E x {GF{A t )-s<H Q <t<H l ) 

(E.18) ft 

= / Eo[F(At-r);t-r<Hi)m XtStG (dr), 

J s 

with the finite, signed measure m XjSi G on being defined by 



(E.19) m XySyG (C) = E x (G;H eC,s<H <H 1 ), C G B{R+). 

Now suppose that k G N, s < t = t\ < ti < ■ ■ ■ < tk, and that hi, hi, ■ ■ ■ , hk are 
bounded measurable functions on [0, 1]. Put 

F(x) = h x {x) {U b t2 _ tl h 2 Ul_ t2 h 3 ■ ■ ■ U b tk _ tk _h k )(x), x E [0, 1]. 

Then equations ( |E.18[ ) and ( |E.19| > entail 

E x (Gh 1 {A tl )---h k (A tk ); s < H < h < t 2 < ■ ■ ■ < t k < H{) 

(E.20) ft 

= J Eo(hi(A tl -r) ■ ■ ■ hk(A tk - r ] t k - r < Hi) m x , s , G (dr). 

Next we make the following choice 

G = g 1 (A Sl )g 2 (A S2 )---g l (A Sl ), 

with I e N, < s% < ■ ■ ■ < si = s, gi, . . . , gi being bounded, measurable functions 
on [0, 1]. Consider the measure m XjS: G in this case, and let C G S(M+). Then 

m x , s ,G(C) = E x ( gi (A Sl ) ■ --giiA^); H G C, s t < Hq < Hi). 

Recall that we assume, as we may, that A is constructed pathwise from the standard 
Brownian motion O, which exclusively has continuous paths. In particular, the paths 
of A and O coincide up to the hitting time Hq (inclusive), so that Hq = Hq, where 
H^ and H® denote the hitting times of c = 0, 1, by A and O respectively. We 
claim that {Hq < H^} = {Hq < H®}. Indeed, this equality is equivalent to 



54 V. KOSTRYKIN, J. POTTHOFF, AND R. SCHRADER 

{H^ < Hq} = {-fff < Hq}, but since A and O coincide on the random interval 
[0, Hq] = [0, Hq] this is obvious. Consequently the formula 

m x>SjG {C) = E x ( 9l (0 Sl ) ■ ■ ■ gi (0 Sl ); H° e C, s < H° < H?) 

holds true. Moreover, starting in 0, A is equivalent to the process A 1 . Hence for- 
mula ( |E.20[ ) reads for the above choice of G 

E x (Gh 1 (A tl )---h k (A tk ); s < H < h < h < • • • < t k < H x ) 

= J" E x ( 9l (O Sl ) ■ ■ -91(01); F °° G dr, s < H°° < H°°) 

x EQ{h x {A\ x _ r ) ■ ■ ■ h k (A\ k _ r ); tk-rKa 1 ), 

and we recall that a 1 denotes the hitting time of 1 by A 1 . From the construction of Y 
we obtain 

E x ( 9l (A Sl ) ■■■gi(A Sl ) 

x h^AtA ■ ■ ■ h k (A tk ); s t < H^ < t x < t 2 < • • • < t k < Hf) 

= E x { 9l (Y Sl ) ■ ■ ■ 9l (Y Sl ) 

x hi(Y tl ) • • • h k (Y tk ); s, < H% < h < t 2 < • • • < t k < Hj) 

= E x { 9l (X Sl )--- gi (X Sl ) 

(E.21) x h x (X tl ) ■ ■ ■ h k (X tk ); siKH* < h < t 2 < ■ ■ ■ < t k < H?) 

and in the last step we used lemma [57T] 

Finally, let f\, f n , n G N, be bounded, measurable functions on [0, 1], and 
assume that < V\ < v% < • • • < v n . Then we write 

E x [fx(A Vl ) ■ ■ ■ f n (A Vn ); v n < H£) 

= E x (h(A Vl ) ■ ■ ■ f n (A Vn ); v n < H* A Hf) 

n 

+ J2Mh( A ^) ■ ■ ■ fn(A Vn ); Vj-! < H^ < Vj , V n < H?). 
3=1 

Since up to the hitting time of 0, A and X are both equivalent to a standard Brownian 
motion, the first term on the right hand side is equal to 

E x (h{X Vl ) ■ ■ ■ f n (X Vk ); v n < H* A iff). 

For each term under the sum we use formula ( |E.21[ ) with appropriate choices of I, k, 
s%, gi, i = 1, . . ., I, tj, hj, j = 1, . . . , k. Therefore we finally obtain 

E x {h(A vl ) ■ ■ ■ f n (A Vn ); v n < Hf) = E x (f x {X vl ) ■ ■ ■ f n (X Vn ); v n < H*), 



and the proof of the first statement of lemma |5.5| is finished. The second statement 
follows immediately, because A is a Markov process, and hence so is A (e.g., J3] 
Chapter III. 3], [6, Chapter X]), and the simple Markov property actually is a property 



of the finite dimensional distributions of a process. The last statement of lemma 5.5 
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is then a consequence of the fact that A and B are equivalent to a standard Brownian 
motion up to their hitting time of 0, 1 respectively. □ 

For t G M + , and k G N, we consider the family Q kt °f events A of the form 
(E.22) A = AinA 2 n{T fe < si}, 

with Ai G J^, 

n 

A 2 =P|{X S , gC,}, 
i=i 

and where n G N, < si < . . . < s n < t, C x , C 2 , . . . C n G B([0, 1] U {A}). Then it 
is a routine exercise to show the following result: 

Lemma E.3. Qk,t is a H-stable generator of the a-algebra D {T k < t}. 

Proof of lemma 5. 7 Take A of the form ( |E.22[ ). Observe that T k < s± entails that T k 
is finite, and thus Xx h G {0, 1}, i.e., Xrp k 7^ A. We have to compute 

E x (f(X s+t ); An{s + t<T k+1 }). 

(Note that A n {s + t < T k+1 } includes the event {T k < t < T k+ i}.) We want to 
apply lemma [ETT] and choose G in equation ( |E.12[ ) as follows: m = n + 1, s m = 

s n +i = s + t, and 

n 

g(X Sl ,. . .,X Sm ) = [] l Cj (X s .)f(X s+t ). 



Moreover, we choose the (p in lemma E.l as the indicator of (s + t, +00]. Then we 
can compute with lemma [E~T] as follows: 



E x (f(X s+t );An{s + t<T k+1 }) 

= E x (<p(T k+1 ) G; {T k < si} n Ai) 

= E x (E x [<p(T k+ i) G; T fc < si I jg) ; A^ 

= E x (E XTh {<p(r + H e )G r ) l [0>Sl] (r) \ r=TkjC={Xt y, Ai). 

As before, we decompose Ai = A1.0 W A^i in J 7 ^ with Ai )C = Ai n {Xx k = c}, 
c = 0, 1. The term involving A^o is equal to 

(E-23) E^E^r + H^Gr) l [0 , si] (r) | r=Tfe ; A 1)0 ). 

We consider the inner expectation for r G [0, s\]: 

n 

Eofair + HjGr) = Eo(j{X a+t - r ) JJ lc,(X s ._ r ); s + 1 - r < #1) , 
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and remark that this expectation only concerns the process X killed when reaching 1. 
Moreover, the event r\j{X s .^ r G Cj} belongs to T^_ r . We use the Markov property 
(cf. lemma 5.5 ) of the killed process relative to this u-algebra, as well as 

{s + t-r < Hi} = {s < Hi o e t - r , t-r < Hi} 

to compute 

E (ip(r + Hi)G r ) 

= E (E (f(X s ) o t _ r - s<H l0 Q t _ r | T t x r ); 

{t-r<H 1 }n{r$ =1 {X aj - r eC j } 

= E (E Xt _ r (f{X s );s<Hi); 

{t-r<Hi}n (n? =1 {X s ,_ r G Cj}) 
= E (j r ; t-r < Hi), 



with 



Jr = E Xt _ r (f(X s ); s < H x ) JJ l C] {X s 

3=1 



(;•) | .,. - E (J r ; t<r + Hi) l [0)Sl ](r) 



\r=T k 



Therefore 

£oMr + ffi)G r )l [0iSl] (r /|r=Tfe 

= E X [J ; T k < si, t < T k+ i | J 7 ^), 

with another application of lemma |E.l| Set J = Jq and insert this into ( E.23 ), then 

^(/(X s+t ); A n {X Tfc = 0} n {s + t < T k+ i}) 

= E X (^E X (J; T k < si, t < Tk+i |-^t[); Ai t o 
= E X (J; {T k <si,t< T k+ i} n Ai, ) 
= E x (E Xt (f(X s ); s < Hi); {t < T fe+1 } n A n {X Tk=0 } 



We add this expression to the corresponding one, where {Xt h = 0} is replaced by 
{XT k = 1}, and find 



E x (f(X s+t ); AD {s + t <T k+ i}) 

= E x (E Xt (f(X s ); 8 < H e ) 



.; {t<Tfc+i}nA). 
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Therefore we have shown (P x -a.s.) 

E x (f(X s+t ); T k <t,s + t< T k+1 | Jf ) 

= E Xt (f(X s ); s < H c ) \ c={ x Tk )* 1 {T k <t<T k+1 } 

= E Xt {f(X s ); H c * <s< H c ) | c=(XTfe) * l {T k <t<T k+1 } 

+ E Xt (/PQ; s < H c A H c ,) | c=(Xt j* 1 {T k <t<T k+l } 
= E Xt (f(X s ); Tx < s < T 2 , X Tl = c) \ c=Xt l {Tk <t<T k+l} 

J k 

+ E Xt (f(X s ); S <T 1 ) l {Tk < t<Tk+l} , 
which is the statement of lemma [5771 □ 



Proof of lemma pT5| We work on the event {X Tk = 0, T k < t < T k +i, X t / A}, 
fcEN . Then 

E x (f(X s+t ); T k <t<T k+1 <s + t\ Jf ) 

= E x (E x {f(X s+t ) | F$ ) ■ T k <t<T k+1 <s + t\ Jf 



(E.24) = E x (<p(T k+l ); T k <t\T t x ), 

with 

<p(r) = Ei(/(X s+t _ r )) l( M+t] (r), r G 
and we have applied lemma [ET] Let A E Q k jt be of the form ( |E.22| ), and as before 



we write Ai o = Ai n {XT k = 0}. We set G = 1a 2 , then with another application of 
lemma [ETT] we get 



E x (<p(T k+1 );An{X Th = 0}) 

= E x (<p(T k+1 ); A lj0 n A 2 n {T fe < Sl }) 

= E x (E x (<p(T k+1 ) G | J-#J; Ai,o D {T k < Sl } 
(E-25) = E x (Eo(<p(u + Fi) G„) l [0 , ai ](«) | u=Tfc ; Ai,o 

Consider the inner expectation value 

n 

E (ip(u + £Ti) G„) l[o, Sl ](«) = E ^(u + Hi) ]J 1 C .(X S ._„)) l [0 , Sl ](u). 

The definition of tp implies that under this expectation we have t — u < Hi, and 
therefore it only concerns the process X killed at 1, together with its lifetime Hi 
(which we recall is a stopping time relative to the natural filtration of X). Therefore 
by lemma 5.5 we can use the Markov property of the killed process to compute this 
expectation. To this end note that si — u < . . . < s n — u < t — u, so that G u is T^_ u 
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measurable. Thus 

E (ip(u + Hi) G u ) l[ , Sl ](u) 

= E a (E (ip{u + H 1 )-t-u<H 1 <s + t-u\ F* u ) G«) l [0 , ai] («) 
= E (E (tp(t + Hi)o 9 t - u ; < Fa o 6 t - u < s \ F?_ u ) 

x G u ; t - u < H^j l[ , Sl ](u) 
(E.26) = E (E Xt _ u (<p(t + H 1 );0<H 1 <s)G u ;t-u<H 1 ^ l [0>si] (u), 

and we made use of the relation 



{t-u < Hi < s+t-u} = {Hi = t-u+Hio6 t - u , < HioOt-u < s, t-u < Hi}. 



We evaluate equality (E.26 1 at u = T k , and apply lemma E.l This gives 



E Q (ip(u + Hi) G u ) l [0 , si ](u) \ u=n 

= E [E Xt _ u (ip{t + Hi);0<Hi<s)G u ;t-u<Hi) l [0 , si] («) \ u=Tk 
E x (E Xt (<p(t + fli); < Hi < s) G; T k < s x , t < T fe+1 | J* ' 



The last expression is inserted into equation (E.25 1. This results in 



E x (ip(T k+1 ); An{I Ti =0}) 

= E x (E x (E Xt (<p(t + Hi); < Hi < s) 



xG;T k <si,t<T k+1 \T£);Ai i0 



E x (E Xt (<p(t + Hi);0<Hi<s); {X Tk = 0,T k <t< T k+ i} n a) . 



With equation ( |E.24| ), this implies that on {X Tk = 0, T k < t < T k+ i, X t / A} we 
have the relation 



(E.27) 



E x (f(X s+t ); T k <t< T k+ i <s + t\ T?) 

= E Xt (ip(t + Hi);0<Hi<s), 
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because A^o G J 7 ^. Using lemma E.l we compute the right hand side of this for- 



mula: 

Extfait + Hx); < Ht < s) 

= E Xt (<p(t + Hi); < H < H x < s) 

+ E Xt (cp(t + H 1 y,0<H 1 <sAH ) 

= E Xt (<p(t + T 2 ); T 2 < s, X Tl = 0) 

+ E Xt (<p(t + T 1 ); Ti < s, X n = 1) 

= ExtfatfiX-r)) | r=T2 ; T 2 < s, X Tl = 0) 

+ Ex^EhtfiXs-r)) \ r=Ti ; Ti < s, X Tl = l) 

= E Xt (E Xi {f(X s ) | 7"*) ; T 2 < s, X Tl = 0) 

+ E Xt (E Xt (f(X s ) | Jg) ; Ti < s, X Tl = l) 

= E Xt {f(X s );T 2 <s,X Tl =0) 

+ E Xt {f(X s );T 1 <s,X Tl =l). 

We have proved that on {Xx k = 0, < t < T^+i, Xt ^ A} the following formula 
holds true i^-a.s.: 

E x (f(X s+t ); T k+ i<s + t\T t x ) 

= E Xt (f(X s ); T 2 < a , X Tl = 0) + E Xt (f(X s ); T x < s, X Tl = l). 



The case {Xx k = 1} is dealt with in the same way, and lemma 5.8 is proved. □ 
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